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ABSTRACT 

The values of the NRCS, o° , obtained by a scatterometer a: 

random variables. They are random variables with the property thi 

the variance is a known function of the expected value (or the me; 
in a loose sense) . The probability density function has featur< 
that allow required probabilities to be obtained from the norm, 
distribution. Models for the expected value obtain it as a functi' 
of the properties of the waves on the ocean and the winds th, 

generated the waves with the goal of obtaining wind speeds a: 
directions. 

A scatterometer presently in the design stage has improved fe 
tures compared to the SASS on SEASAT in that the backscatter valu 
are closely located on a well defined grid. Also the requir 
parameters are known from the design features of the instrument. 

Point estimates of the expected value, c° , are found fr' 

various statistics given the parameters that define the probabili 
density function for each value. Random intervals are derived wi 
a preassigned probability of containing the value, o°. 
statistical test to determine whether or not successive values 
o° are truly independent is derived. 

Ways to correct for certain kinds of errors in a model sole 
by means of backscatter data are derived. Ways to use convention 
measurements from the National Data Buoy network to interpr 
backscatter models more thoroughly and to either improve them 
validate them are found. 


vi 


The maximum likelihood estimates for wind speed and direction 
are found, given a model for backscatter as a function of the 
properties of the waves on the ocean. These estimates are biased 
as a result of the terms in the equation that involve natural 
logarithms, and calculations of the point estimates of the maximum 
i likelihood values are used to show that the contributions of the 
logarithmic terms are negligible and that the terms can be omitted. 
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1 . INTRODUCTION 


1.1 Historical review . The term "scatterometer" was coined 
by R. K. Moore (Personal Communication, 1986) and used in the de- 
scription of a radar-radiometer by Moore and Ulaby (1969) . A 
scatterometer differs from other radars in that estimates both of 
the received signal power plus the receiver noise power and of 
the receiver noise power alone are obtained. The latter is then 
subtracted from the former to obtain an estimate of the received 
signal power alone. Both of these randomly varying quantities 
are measurable and of course, by definition, positive. However, 
if the randomly varying received signal power is very small com- 
pared to the randomly varying receiver noise power, subtracting 
the randomly varying noise power from the randomly varying signal 
plus the noise power can produce a negative value for the ran- 
domly varying estimate of the received power. 

Scatterometer data are used to attempt to determine the 
speeds and directions of the winds over the ocean by means of 
backscatter estimates obtained from a scatterometer on a polar 
orbiting spacecraft at an altitude of about 800 km. The first 
attempt to do this globally was by means of the SASS on SEASAT. 
At first it was believed that the wind speeds and directions ob- 
tained from the backscatter estimates were quite good as in Duffy 
and Atlas (1986), Pierson, et al. (1984), and Woiceshyn, et al. 
(1985). Further studies have revealed various kinds of 
systematic errors in the wind recoveries, most noteably as de- 
scribed by Woiceshyn, et al. (1986) . 
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1.2 Design features . Schematic versions of a new scanning 


pattern are shown in Figs. 1 and 2. If to scale, each square 
would represent an area of the ocean surface that is 25 by 25 km. 
Two possibilities for obtaining backscatter data for each of these 
areas are under consideration. 

The scanning pattern as in Fig. 1 consists of 25 by 25 km 

squares oriented both parallel to the subsatellite track at 

nearly a constant, but different, incidence angle for each of 

the four measurements over considerable distances and normal to 

o o 

the track for incidence angles that vary from about 15 to 59 . 
The along track variation of the parameters associated with each 
backscatter value is small over distances of several hundred km, 
but the cross track variation is rather large. There are 24 
areas, each 25 by 25 km, for a total distance of 600 km normal to 
the subsatellite track on each side of the track. The incidence 
angle, including the beam at 115°, varies from about 15.6° to 
59.1°. 


One possibility, as shown in Fig. 2, is to have one ver- 
tically polarized and one horizontally polarized value for the 
backscatter estimates obtained from forward looking antennas at 
45° to the subsatellite track to the right side. A third 

estimate would be obtained by means of an antenna at 115° to the 
track, quite possibly for horizontal polarization. The fourth 
for vertical polarization would be at 135° to the track. 

The second possibility is a vertically polarized estimate at 
45°, both vertically and horizontally polarized estimates at 115 
and a vertically polarized estimate at 135°. The first of these 


i 


2 



>0 



Fig. 1 Schematic Drawing of Large Scale Grid of 25 by 25 Km Areas 
The Chessboard Pattern in the Expanded View Actually Fills 
the Entire Area. All Squares Provide Data. 


x-x 0 




Fig. 2 Two Possible Combinations of Four Backscatter Estimates 
for the Right Side Looking Forward. 
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two possibilities seems to offer the surest opportunity to 
recover unique values for the wind speed and direction because 
of the apparent marked difference between upwind backscatter 
estimates and downwind backscatter estimates when vertically and 
horizontally polarized values are compared for high winds and 
high incidence angles. The first of these two possibilities will 
be used in a subsequent section. The changes required for the 
second possibility are easy to make. 

Fig. 2 illustrates a feature of the actual data in that the 
pointing directions of the radar beams are not exactly 70° and 
90° behind the lead beam. This produces minor shifts in the 
various inverse curves that result and complicates the recovery 
of the winds from a computational point of view. This aspect of 
the problem will be neglected in the analyses and discussions to 
follow. 

This new design as in Chi, et al. (1986) and Chi, et al. 
(1987) and in terms of various other sources of variability for 
the estimates, such as attitude errors, orbital errors, and 
squint, is a vast improvement over the SASS. The major dif- 
ficulty will be to obtain a correct theory, or model, to relate 
the backscatter estimates to the waves on the ocean and then the 
winds that generated the waves. 

1.3 Probability concepts . Both the backscattered signal 
from the waves on the ocean and the receiver noise are in- 
coherent, which means that the envelopes of the time varying 
signal plus the noise as, say, voltages can be described instan- 
taneously by a Rayleigh probability density and that the power 
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has an exponential distribution. Successive values from these 
envelopes are correlated for short time intervals. 

The Doppler shift as a function of the position of a given 
small area of the ocean and of the velocity of the spacecraft al- 
lows successive small areas as in Figs. 1 and 2 to be sampled and 
estimates of the backscatter to be obtained for each one. The 
new feature of the present design is that digital data processing 
and filtering will make the precise alignment of the areas il- 
lustrated in these two figures possible. 


Probability theory and statistical methods are involved in 
the design of a scatterometer as described by Chi, et al. (1986) 
and Chi, et al. (1987) . The design involves combining a large 
number of guasi independent estimates of the random variables ob- 
tained by the instrument so as to obtain a much more stable, less 
randomly varying, estimate of the backscatter. 


Although the theory of sampling from a normal distribution 

with an unknown expected value for the first moment about the 

. . . 2 2 
origin, conventionally, y, and an unknown variance, £ (x-vO =o 

(note this is not o°) , does not quite apply, a scatterometer es- 
sentially obtains a large number of estimates of a random 
variable, say, X Jf X 2 , X^, .... X n and finds the average as in 



( 1 ) 


It then, in a way, follows that the expected value of X is 
given by 

e (X) = y (2) 
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and that the variance of x is 


j e(X - y) 2 = a 2 / N (3) 

The single number that is obtained from one of the four dif- 

i 

ferent estimates of the backscatter from the waves in a given 25 
by 25 km area on the ocean is thus analogous to the X in Eqn. 
j (1), and the variance of the estimate (i.e. X ) is analogous to 
the expected value given in Eqn. (3) . The analogy is not exact. 

; Therein lies both the value of a scatterometer and the dif- 
ficulty in interpreting the estimates of the backscatter. One 
purpose of this paper is to explain this difficulty and to 
provide a correct theory for the interpretation of the 
backscatter estimates. 


1 
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2 . DEFINITIONS 


2.1 The quantities involved . For the analysis to follow, 
certain definitions are needed. These are: 

o° is a random variable. 

o o 

a M is the expected value of o . 

VAR (a 0 ) = a (a. 0 .) 2 +6 a° + y is the variance of o° . 

SD(o°) = (VAR (o° )) 2 is the standard deviation of o°. 

Since the variance is a known function of o° > VAR(o°) and 

SD ( o°) will be replaced by VAR(a°) and SD (c°) throughout. 

a, 6 and y are known design parameters for the scatterometer 

for a known 25 by 25 km cell and polarization. 

o 

In these definitions, a is the value of the normalized 
radar backscattering cross section obtained by a scatterometer, 
and is its expected value. The purpose of the definitions is 

to tie together concepts used in the theory of probability and 
statistics with those used in scatterometery and to avoid confu- 
sion between the notations used in the two different subject 
areas. 


2.2 The probability density function . With these 
definitions, the probability density function for o° is given by 


f (0°) = (2u)' }2 (SD(a°))' 1 exp (- j(a° - c°) 2 /VAR(a°)) (4) 
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f 

as first obtained by Fischer (1972). It is not usually necessary 
to designate the polarization of the backscatter in the following 
analysis so that the appropriate subscripts will be omitted until 

they are needed. There is only one unknown parameter in Eqn. 

o o 

(4). It is o M ,and if it is known, both the expected value of o 

and the variance of o° are known. The range of definition for o° 
is 

- CO < 0° < 00 

( 5 ) 

but the range of definition of cr° is 

0 < o° < 00 

— M 

because the theory for backscatter from any target does not admit 
- negative values. It is conceivable, however, that there could be 
no backscatter at all from certain kinds of targets, or that the 
value of a°j could be so low that it could be virtually in- 
distinguishable from zero. Eqn. (4) still has a meaning if 0 ° 
is zero. For the derivation of Eqn. (4) , it is assumed that both 
the expected value of the receiver noise and the variance of the 
[ noise are constant. The backscatter estimates recovered by the 
j SEASAT SASS for some incidence angles ranged over 40 dB. For 
some theoretical purpose it may be necessary to show results in 
dB if o° is positive. For the practical uses of the data, there 
does not seem to be any reason to use dB for the recovery of 
winds. 

i 
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3. 


THE MEANING OF a° 

M 


3.1 Backscatter theories . Theories for the calculation of 
radar backscatter from the ocean surface attempt to predict the 
value of from either the properties of the waves or the study 
of backscatter data. The most recent models are those of Wentz, 
et al. (1984) and Wentz, et al. (1986), Iwata (1985), Durden and 
Vesecky (1985), Plant (1986) and Donelan and Pierson (1987). 
These models do not treat separately some properties of the waves 
known to produce backscatter such as breaking waves (Banner and 
Fooks, 1985) and wedges (Lyzenka, et al. 1983, Kwoh and Lake, 
1984). In general, o° is some function of the properties of the 
waves and thus indirectly some function of the winds (both past 
and present over the entire ocean) that generated the waves. 
More specifically, it should be possible to represent o° as 
dominantly a function of the wind speed and the direction of the 
wind (relative to the pointing direction of the radar beam) for 
the area being analyzed and to account for such effects as fetch, 
duration, swell, breakers, wedges, and water temperature in the 
process of estimating the winds. To do this requires the com- 
bined efforts of oceanographers, meteorologists and radar scien- 
tists. Numerous estimates of backscatter obtained by measure- 
ments from radars on platforms in oceans or lakes, on aircraft 

and on SEASAT have yielded data that can be compared with models 

_ o 

for a .. . 

M 

The last four models referenced above (there may be others) 
are based on the two scale theory reviewed by Valenzuela (1978) . 
The separation of scales wavenumber is a matter of considerable 
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debate. Backscatter theories that claim to avoid the separation 
of scales assumption have both proponents and opponents as in 
Bahar (1987) and Guissard (1987) and the references contained 
therein. 

One purpose of this paper is to find bounds on so that 
various models can be tested operationally. The value of a° for 
a particular sample value of a° from Eqn. (4) can be assumed to 
be known theoretically based on some model for backscatter from 
the waves on the ocean. For purposes of analysis in the material 
on maximum likelihood estimates and the use of various models it 
will be assumed that, with polarization included, 


' °W<"' x - e) (7) 

°HH ' °HH (0 - X> 6) <8> 

where 0 the incidence angle, is known. 

3.2 The specification of the wind . It is not necessary to 
specify the height at which the mean effective windspeed 0 for a 
neutrally stratified atmosphere is to be measured. Some inve- 
stigators use u* instead of 0. If the logarithmic wind profile 
for a neutrally stratified atmosphere is used and if the friction 
velocity is determined by one of the usual methods to be 


u. 


(O(io)y 


c dn W 10 )) 


(9) 
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then u 4 and z q are both functions of 0 (10). The wind at 10 m, and 
the wind at any height, or u* are all equivalent. In the Monin- 
Obukhov theory, all wind profiles are logarithmic close to the 
sea surface. There are differences as far as some backscatter 
theories are concerned, and there are many different versions for 
the right hand side of Eqn. (9) . Also for other than neutral at- 
mospheric stability, Monin-Obukhov theory is needed to interpret 
comparions of conventionally measured winds with backscatter 
data. The air-sea temperature difference is required from some 
other source to compute the actual wind at greater heights from 
the effective neutral wind. Hasse (1986) has discussed the 
various concepts that are involved. Donelan (1982) has provided 
an analysis to support the contention that the right hand side of 
Eqn. (9) is also a function of the waves on the ocean. The 
values of (7) and (8) , or improved models for o°, vary over many 
orders of magnitude in such a way that it should be possible to 
recover estimates of 0 and X from the backscatter estimates. 
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4. PROPERTIES OF EQUATION 4 
4.1 Monte Carlo methods . The transformation 

t = (a° - a°)/(a(a°) 2 + $ 0° + Y ) h {10) 

changes Eqn. (4) to Eqn. (9) . 

f(t) = exp(-t 2 /2) (11) 

The random variable, t, is a sample value from a normal 
probability density function with an expected value of zero and a 
standard deviation of one. From (10) it follows that 0° can be 
represented by Eqn. (12) . 

0° = 0° + t • SD(0°) (12) 

The sampling variability of 0 ° in Eqn. (4) is entirely equivalent 
to the values obtained fron Eqn. (12). If some model for o° is 
available, even if it is not fully correct, then can be used 
to generate random values of o°. These in turn can be used to 
obtain winds recovered from simulated o° estimates and test the 
algorithms designed to do this. The theoretical values of 

need not necessarily be correct for this to be done. The model 

values, 0 ^, need not necessarily correspond to the best model 
values that would recover more nearly correct winds. 
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Eqn. (12) can 


4.2 Kp The normalized standard deviation , 
be used to obtain K p as follows: 


o° = o° (1 + t • SD (a°)/a°) 


(13) 


= °M ^ 


+ t 


K p )) 


(14) 


As a° approaches 0, K p approaches infinity. As o° approaches in- 
finity, K p approaches a 2 . 

4.3 The normalized expected value . From Eqn. (10) if o° is 
known, a random variable can be defined by 

Y = C°/SD(0°) 

and its expected value is 

e 00 = m = °°/SD(o°) 

The probability density function is then 

f(Y) = (2tt) exp (- %(Y - M) 2 ) (17) 

The variance of Y is one. Eqn. (17) has the same shape as the 
graph of the unit variance zero mean normal distribution, but it 
is centered at the value, M. For 0 ° = 0, M = 0 . As 0 ° becomes 

M M 

. -h 

large, M asymptotically approaches a . Also M is the reciprocal 
of K p defined as in Eqn. (18). 

M = 1/Kp (18) 


(15) 


(16) 
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| Given a°, a, 3 and y the probability that Y and hence c° 

| will be negative can quickly be guessed, or computed, since M 
[ represents the number of standard deviations that the expected 

value lies away from zero. For large values of a° and for one of 

the areas that would be sampled, a possible value of a is 0.02 so 

that a 2 equals 0.14 or 7.07. The integral of the unit 

i 

variance zero mean normal probability density function from-°° to 
-7 is a very small number. Thus a negative estimate of o° would 
be highly improbable. For a° equal to zero, M is zero, and the 
probability of a negative value of Y or o° is 0.50. The standard 
deviation of a is then y 2 . 

t 

Although Eqn. (4) has a range for o° over all possible posi- 
I tive and negative values and for o° from zero to infinity, 
probability models are both unrealistic and unverifiable over 
their entire domain of definition. The actual instrument sets 
physical limits that constrain the random variables and unknown 
parameters in ways that are features of its design. A value of M 
equal to 4 or more is sufficient to ensure that the chance that 
o° will be negative is negligible. 

* 
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5. 


POINT ESTIMATES OF POPULATION PARAMETERS AND RANDOM 
INTERVALS FOR POPULATION PARAMETERS 

5.1 The concepts . Given an independent random sample of 
values from a population with a probability density function with 
a known form, but with one or more unknown parameters that define 
it more exactly, the classical statistical problem is to estimate 
the unknown parameters. An independent random sample from a 
normal distribution with an unknown expected value and an unknown 
variance provides the data from which various statistics can be 
computed so as to obtain point estimates of the two parameters 
and so as to construct intervals whose end points have a pre- 
assigned probability of enclosing the values of each of the two 
population parameters. The procedures are a part of the clas- 
sical literature on the subject and the random intervals are con- 
structed from statistics that involve the "student t" distribu- 
tion and the Chi Square distribution. All terms and the methods 
that are used are carefully defined in, for example, Mood, et al. 
(1974) , except that the intricate notation of this reference has 
been simplified in this paper. 

For scatterometers in general and for any new design, a 
statistician is thus faced with the task of stating as much as 
possible about the unknown population parameter, o° , given at 
first a single value, i.e., a sample of size one, for the random 
variable, o°. Clearly the standard methods applied to a sample 
from a normally distributed population do not apply. As de- 
scribed above, the data from the scatterometer have already been 
processed extensively to the extent that the random variable, 0 ° / 
is like the x in Eqn. (1) and the variance is like Eqn. (3) for 
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quite large n with the major difference being that the variance 
is a function of the first moment about the origin. 

5.2 Probabilities . For purposes of further analysis the 
random variable t in Eqns. (10), (11), and (12) will be, at 

times, defined with subscripts so as to differentiate different 
values. Also at times the absolute value of these subscripted 
values will be assigned a sign to designate whether the sub- 
scripted value is positive or negative. 

The integral of Eqn. (11) is often given in tabular form 
(pre-computer) from minus infinity to some positive value of t, 
say, +|t 2 | , as in 

i r |t2 ' 

$(+|t 2 |) = (2 tt) exp(-t 2 /2)dt (19) 



I so that 

f+|t 2 l 

$(+|t 2 | - $(-|t 1 |) = (2 tt)"^ exp(-t 2 /2)dt 

Hhl 

With these definitions, it follows that 


( 21 ) 
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P(o° - 1 1 1 1 SD(o°) < O°<0° + |t 2 |SD(o°)) = $(+|t 2 |) - $( -Itjl) (22) 

Where F ( • ) means that the probability that o° will fall in the 
interval defined by - 1 1 | , o° , SD(o^) and +|t 2 | is given by the 
right hand side of the equation. 

5.3 Point estimates . Given a sample from a known 
probability density function, hereafter PDF, with unknown 
moments, the first moment about the origin can be estimated by 
the sample mean. The average of one number is that number itself 
so that 

o o 

- ° (23) 

is an estimate of o°. This estimate is not a particularly good 

M 

one since, from Eqn. (22), it scatters by a considerable amount 
about the expected value of o° for, say |tj| = |t 2 | = 1 . This 
point estimate was substituted into the equation for the standard 
deviation to estimate the scatter of o° for reduction of the 
SEASAT SASS data. 

An alternative point estimate is the maximum likelihood 
estimate of a° . For a given sample o° is known, and the maximum 
likelihood estimate of o° is that value of o° that maximizes 
f (a 0 ) for the sample value of o°. Omitting the ( 2 tt ) -35 , multiplying 
by, or adding, constants does not change the location of the max- 
imum. Also since the natural logarithm of f ( o°) is more con- 
venient to work with, the result can be put in the form of Eqn. 
(24)* 
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I 


- 0 , 


MLE (o.J = MAX 
M 


-h(- 


M 




-h 2n (- 


VAR(a°) 


a(o°) 2 + B o° + y 


(24) 


as 0° is varied. The constant in the denominator of the 

logarithmic term causes the function to be zero when o° = o° . 
The derivative of Eqn. (24) with respect to o° when set to zero 
gives the MLE if that part of the PDF is considered for which 

o° f 0 . The result is a cubic equation. A simple search by a 

M 

o 

computer subroutine can generally locate MLE( ° M ) . It is 

slightly less than . Values of the MLE will be given later 

based on newly available design parameters. 


If both mesoscale variability from one 25 km resolution 
area to the next (Pierson, 1983a) and the sypnotic scale 
gradients can be neglected in the sense that both are small com- 
pared to the sampling variability of estimates of the backscatter 
(Pierson, et al. 1984) when various numbers of 25 km resolution 
cells are combined, it should be possible to average the 
estimates of o° to obtain more stable estimates. 


o 

Consider for example, n values of 0 for a 

. o o o 

the subsatellite track as in a, , o 0 , ... o . 

f 1 *• Ti 

t would be 25 km wide by 25 m n km long. Then 

t 


-o 

a 


i n 

1 V o 

"i °i 


1 , o 

“ ( n °w 

n M 


* l 


row parallel to 
The area sampled 


(25) 


from Eqn. (12) . The expected values of o° and (o° - o°) 2 are 
then 
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p r O 

e (a ) = o M 


( 26 ) 


and - a°) 2 = (i| t. SD(o°)) 2 

= ijn VAR(a°) 
n 

VAR(a°) (27) 

n 

2 

since &(t i ) = l and e (t , t ) = 0 (p^q) . Effectively, averaging n 

values of o° along a row reduces a, 6 and y by n” 1 . 

Similarly, areas that are n long by n wide could be repre- 
sented by 



o° . + t. . SD(o° .) 
Mi} ij Miy 


(28) 


The average would be approximately 


(29) 



so that 

£(a°) = a° = effective value of the (30) 

where 5^ could represent the effective model value at the center 
of the larger area. 

The variance would be 

e < 5 ° ' ^ = 7 I f VAR(5 M ) (31) 


20 


! t 2 ♦ 

The a s, 3 s, and y s for the n areas could be summed to find an 

-2 

effective value for them. The variance would be reduced by n 
Further analysis is needed to decide on some of these steps. 

5.4 Interval estimates . The concept of an interval 

estimate is described in Mood, et al. (1974) . Fig. 3 is an il- 

. o 

lustration of this concept. The horizontal axis is and the 

vertical axis is o° for the ranges defined by ( 5 ) and ( 6 ). For, 

say, |t | = |t | t the two continuous curves in Eqn. (22) can be 

o o 

graphed. The range of o for a particular <? M is shown by the ver- 
tical line, and o° is a function of a.°. . 

The probability that a° will equal a value somewhere on this 
vertical line is given by the right hand side of Eqn. (22) . For 
any of these possible random values of o° that lie on the example 
vertical line, a horizontal line can be drawn between the upper 
and lower curves. Any one of these lines will intersect the ver- 
tical line for the example value of a° . Two examples, out of an 
infinite number are shown by a continuous line and a dashed line. 
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Fig. 


0 ° - <$( 0 °. | tl |) 


(33b) 


, which is equivalent to interchanging the independent and depen- 
dent variables, then 


PC 0 °( 0 °. |t 2 |) < o° < o°(a°, |tj|)) - *(*|t 2 |) - ♦(-lt 1 |) < 34 > 

if a°(o°,|t 2 |) is greater than zero. 

Eqn. (34) describes an interval for this case with lower and 

o 

upper bounds that are a function of the random variable, o , 
i.e., the measurement made by a scatterometer, which is known. 
The value of o° is not known yet. The interval is consequently a 
statistic. It can be computed, given a, 6 and y for each sample 
value. The parameter, a° , has the above probability of being con- 
tained in this random interval. 


If c° > + |t 2 |y' s as shown by a dash-dot line three additional 
probabilities can be found. The value 


- |t 3 | = o°/y h 

can be computed. Then 


(35) 


P(0° f 0) = l - *c+ |t 3 |) 

(36) 

p (° < o° < a°(o, jt 2 |)) = *(+|t 3 |) - $(|t 2 |) 

(37) 


and 
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I 


P(o°(a°, ItJ) < a°) = *(- |t 1 |) 


(38) 


The sum of the probabilities in (34), (36), (37), and (38) 

o i- 

is one as a check. The equations divide the area above a =+|t 2 |y 5 
in the first quadrant into the line a° = 0 and region I, region II 
and region III as shown. 


0 

If a falls between + |t 2 |y 2 and 
and V in Fig. 1 the value of, say, 



as in Regions IV 


t 4 = 


o°/y h 


can be found even if a° is negative, 
defined as follows: 


(39) 

Three probabilities can be 


P(o° = 0) = 1 - *(t 4 ) 

(and, in particular, if 
is 0.5) 


(40) 

o° is zero, the probability that o° f 0 


p(0 •= °m < °m l<3 °- i'l 1 ” ■ *^4)- *(- lyl) 

P(o°(o°, |tj|) < 0°) = H- |tj|) 

Again the probabilities add to one. 

o i- 

If o is less than - 1 1 1 1 y 2 as in Region VI of Fig. 


(41) 

(42) 


1, the 


quantity 

O » ^5 

t 5 = a /y 2 (43) 

can be found. The two probabilities are 
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! P(c° = 0) = 1 - $(t ) (44) 

t M •> 

i 

• and 

i 

P(0 < 0°) = $(t c ) (45) 

M z> 

It is important to be able to interpret negative values of 
the random variable, o°. The triangular area in Region V of Fig. 
1 shows that there is a non-negligible probability for a random 
interval that will contain a positive value of even if o° is 
negative. 

Two typical choices for |t 2 | = 1 1 1 1 are 1 and 1.645. For 
the value, 1, the right hand side of (34) equals 0.6826. For the 
value, 1.645, it equals 0.90. There is no particular reason why 
four different curves could not be graphed on Fig. 3 to delineate 
these other probabilities. 


6 


SIMPLIFIED RESULTS FOR LARGE VALUES OF a° 

M 


The standard deviation of a° can be re-written as Eqn. (46) 

where 


SD C<$ 




B_ 

2a 


(46) 


For large values of o° the first term in (46) dominates the 

M 

second, and the second term can be neglected to obtain an approx- 
imate result as in Eqn. (47). 


P (°M " i t i i a ^( a M + e/2 < ° C 


" + 


t 2 \ a2 (°M + e/2 °0) 


= h+ ]t 2 o - n- |t 1 i> 

The inequality, 


(47) 


o°{l - ItJ at*) - |t 1 | 6/2 h < o° 

is equivalent to 

a 0 + | t 1 | B/2a^ 

_o , l 

a.. < n 


(48) 


(49) 


with a corresponding result for the other one, so that 


( a 0 - | t. | B/2a^ a 0 + | t | B/2 a^ 

p f . < a < - j 

^ 1 + |t 2 l a^ 1 - 1 t 2 | 


(50) 


-*(+|t 2 |) -♦(-Itjl) 


26 


I 


\ 

( Eqn. (50) defines a random interval, by two statistics such 

i 

j that the probability that this random interval contains the 

| population parameter, is given by the right hand side of the 

1 equation. By assumption, the probability that o° is zero is 

I M 

negligible. Eqn. (36) is then equal to zero. Eqn. (37) equals 

1 - $(+ jt 2 () and so on. 
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7. 


THE SOLUTION FOR ALL VALUES OF 0° 


For the general solution, one needs to replace Eqn. (22) by 
inequalities as functions of o that bound c? M as in Eqns. (32a) 
to (34) . The algebra is best carried out by setting the ine- 
quality sign to an equality sign, squaring an appropriate expres- 

. . o 

sion, solving the quadratic equation that results for o M as a 

function of a 0 , a, B, y and |t.J or |t 2 | and picking the correct 

sign for the solution. The two results are 


o 0 +|t 2 ^B/ 2 -((o 0 +|t 2 |^ 8 / 2 ) / -(l-|t 2 ^a)((a 0 ) 2 -|t 2 | z y)) 

1 - |t 2 | 2 a 


aSco°.|t: 2 U - 


(51) 


and 


n _ o°+|t J 2 B/2 + ((o°+|t | 2 8/2) 2 -(l- 

o°(o 0 ,|tJ) = 


|t 1 | 2 a)((o 0 ) 2 -|t 1 | 2 Y)) Js 


1 - 


Itj'a 


(52) 


so that the terms in Eqn. (34) become known. Again the random 
values for the bounds for the inequalities are statistics that 
can be computed from the sample value of o° . 

The terms l - 1 1 1 1 a and 1 - |t 2 | a would cause difficulties 
if either one were to become negative. The mathematical model 
for the PDF then fails to be realistic. As discussed above the 
value of a is so small that this does not affect the usefulness 
of the results. 
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For Eqn. (51), if o° = |t 2 |'y' 2 , a° (a°|t 2 |) is zero and is 

undefined if o° < |t 2 |y^ . For Eqn. (52) if o° = - |t 1 |y^ the 

numerator becomes 



(53) 


so that o°(o° 1 1 1 j ) is undefined for o° < - 1 | y' 2 . Fig. 3 shows 
these points and illustrates the requirement that no value of a° 
can be negative. 

These results have many applications to the interpretation 
of backscatter estimates from scatterometers, as will be il- 
lustrated by specific examples in following sections. An example 
of their implications is to consider some typical numerical 
values. Suppose that y = 10 6 and two estimates of o° have been 
obtained: one is equal to 10 5 and the other is equal to -10 5 . 
From Eqn. (39) for the first value of a° , t 4 = o.Ol, and for the 
second, it equals -0.01. 

For a° equal to 10 3 , the following probability statements 
concerning random intervals can be made if 1 1 1 1 = |t 2 | = 1 

P(a° = 0 ) = 0.496 (54a) 

P(0 < °M < °M (a °’ = °* 3453 < 54b > 

P(a° (a°, | tjl) < o°) = 0.1587 ( 54c ) 

The first probability statement would need to be interpreted here 
and when first described above quite carefully. This interpreta- 
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tion would be: Given a large number of sample values of o° equal 
to 10" 5 such that y equaled 10” 6 and 1 1 1 1 equalled 1, drawn from 
populations with unknown values of cr° some of which were iden- 
tically zero, then the probability that 0 ° was drawn from a 
population such that 0° was in fact zero is 0.496. 

If on the other hand, the sample value of o° is -10 5 the 
probabilities become 


P(0° = 0) 

= 0.5040 

(55a) 

P(0 < a° < 

a°(a°, Itjl)) = 0.3373 

(55b) 

P (°M ( °°’ 

ItJ) < 0°) = 0.1587 

(55c) 


where, for known values of a, 3 and y the numerical values of the 
statistics inside the probability statements have changed. It 
would obviously be difficult to distinguish possible values of 
o° one from the other, if successive sample values for o° were 
10" 5 and -10' 5 . 

If, to continue, 0 ° were to equal -1.645*10 3 (the corres- 
ponding positive value would be about -29.8 db) , and j 1 1 1 were 
1.645, then the probability that o° would be greater than zero 
would be 0.05. 

From Eqns. (25) to (31) an average of four values of o° ef- 
fectively reduces the values of a , 3 and y by a factor of four 
and the corresponding standard deviations by a factor of 2. The 
vertical line for the example value of o° would be half as long 
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in Fig. 1, and Region II would be half as wide. The random inter- 
vals obtained above would have their lengths changed. Those cor- 
responding to Egn. (34) would be about half as long for the same 
probability on the RHS. For averages of 16, if possible, the in- 
tervals corresponding to Eqn. 34 would be about one fourth as 
long. 
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8. EXAMPLES BASED ON A NEW SCATTEROMETER DESIGN 


8.1 Theoretical Considerations . There are two sources of 
variability in the backscatter values that could be obtained by a 
scatterometer . These are (1) the variability in the estimates of 
the received power as a result of the combined effects of self- 
noise in the receiver and Raleigh-fading for the backscattered 
power and (2) the uncertainties in the calculation of the rest of 
the terms in the radar equation that convert the estimate of the 
backscattered power to the estimate of the normalized backscat- 
tering cross section, o°. 

The theory for the variability of the estimates of the 
received power is well in hand, but the calculation of the ef- 
fects of the rest of the terms in the radar equation poses some 
important questions. 

The radar equation can be simplified for the purposes of 
this analysis to the form given by Eqn. (56) 

a° . P R R 

In Eqn. (56), the random variable, P R , is well 
from a probabilistic and statistical point of view, 
variable, R, is not. This term, R, depends on the calculation of 
the spacecraft's position, velocity and attitude at each time 
that an estimate of a 0 is obtained. The spacecraft may be a 
little bit ahead, or behind, a little bit to the right, or left, 
and a little bit higher, or lower, than the calculated position 
of its center of mass. It may be traveling a little bit faster 


(56) 

understood 
The random 
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or a little bit slower than calculated. Also the calculated 
angular displacements in roll, pitch and yaw may be a little bit 
high or a little bit low. 

All of these uncertainties enter into the calculation of R, 
and, according to present analysis methods, result in additional 
variability in the estimates of o°. For analysis purposes, as- 
sume that (56) can be represented by two parts as in Egns. (57) 
and (58). 

P R - P M * V D < P M> (57) 

R “ «m * V Dt V < 58 > 

In Eqns. (57) and (58), t p is known to be a normally 
distributed random variable with a zero mean and unit standard 
deviation. The nature of the estimation process for the received 
power assures that each and every value of P R , and hence t p , is 
totally uncorrelated with all other values. It is generally as- 
sumed that t R has this same property in that it would also be a 
normally distributed random variable with a zero mean and unit 
variance. 

However, t R need not necessarily have the property that suc- 
cessive values of t are independent in the probability sense. 

K 

As reviewed by Douglas, et al. (1987) who describe difficulties 
in calculating the orbit of a satellite, the residual errors in 
the calculation of R could be systematic, but unknown, and slowly 
varying along the orbit. About 2.5 seconds is required to obtain 
the eight measurements to be grouped eventually to form the four 
backscatter values estimated for the 25 by 25 km areas on 
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each side of the subsatellite track. There is one possibility 
that t R may be substantially constant over the time interval from 
T 1 to T 3 2 so the unknown errors may be highly cor- 
related from one estimate to the next. A partial correlation is 
also possible. The present probability model assumes that succes- 
sive values of t are independent. 

Eqns. (56), (57) and (58) combine to yield Egn. (59), 


0 ■ Wm * t p SD f p M ) R M * W f K ♦ Vr sd(p m> sd ‘V 


(59) 


which can be rewritten as Eqn. (60) if the last term is omitted. 


°° = °M + + *r sd (°S) 


(60) 


From (60) , the following can be derived. 


e (o°) = 


(61) 


C(c° - 0 °) 2 = e (t p SD(o°) + t R SD(a °)) 2 

= VAR(a°) + VAR(o°) 

P K 


(62) 
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8.2 A noise only model . Table 1 has been summarized for 10 
selected incidence angles about 5° apart across the swath. The 
table contains the values of a, 8 and y as defined in Section 2.1 
that would be appropriate if R in Eqn. (56) were a constant. 
Other quantities of interest are also tabulated. For example, if 
|t 2 | = 1 1 ^ | = 1 with 8/2 a 2 neglected and for high values of the 
backscatter, the random interval defined by Eqn. (50) can be 
found approximately in decibels in terms of the first two columns 
after the tabulated values of a . Receiver noise and Rayleigh 
fading alone make the random interval based on the estimate, cr°, 
at most + 0.25 dB and as small as + 0.12 dB. In what follows the 
square root of y will also be of interest. The last three 
columns provide information on the values of o° such that the 
asymptotic behavior can be used to sufficient accuracy. 

8.3 Position, velocity and attitude errors . The effects of 
the errors in the calculation of R when combined with the types 
of variability represented by Table 1 are shown in Table 2. If 
a 1 , Bj and Yj are from Table 1 and a 2 , B 2 and Y 2 are from Table 
2, then the second term in Eqn. (62) would be given by Eqn. (63), 
which follows from the properties of normally distributed random 
variables. 

VAR(o°) = (a 2 - otjHo®) 2 + (B 2 - 3^0° + (y 2 - Yj) (63) 

A comparison of Tables 2 and 1, shows that the largest changes 
are in the values of a 2 compared to . The ratios o^/aj for the 
incidence angles in the table are given by the values in paren- 
theses following the incidence angle as: 15.6° (10.83) , 
19.9 °(22.88) , 24.7° (22.55) , 31 . 4 0 ( 19 . 42 ) , 35 . 4° (13 . 95) , 
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TABLE la SELECTED SCATTEROMETER DESIGN VALUES FOR COMMUNICATION NOISE AND RECEIVER 
NOISE (ACTUAL VALUES HAVE BEEN MULTIPLIED BY THE POWER OF 10 SHOWN). 
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TABLE lb, SELECTED SCATTEROMETER DESIGN VALUES FOR COMMUNICATION NOISE AN 
RECEIVER NOISE (ACTUAL VALUES HAVE BEEN MULTIPLIED BY THE POWER 
OF 10 SHOWN) . 
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TABLE 2a SELECTED SCATTEROMETER DESIGN VALUES FOR THE COMBINED EFFECTS OF COMMUNICATION 
NOISE, RECEIVER NOISE, POSITION, VELOCITY AND ATTITUDE ERRORS. (ACTUAL VALUES 
HAVE BEEN MULTIPLIED BY 10 TO THE POWER SHOWN) . 
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TABLE 2b, SELECTED SCATTEROMETER DESIGN VALUES FOR THE COMBINED EFFECTS OF 
COMMUNICATION NOISE, RECEIVER NOISE, POSITION, VELOCITY AND 
ATTITUDE ERRORS. (ACTUAL VALUES HAVE BEEN MULTIPLIED BY 10 TO TH 
POWER SHOWN) . 
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50.5° (11.52) , 


and 


40.0 °(12.46) , 46° (10 . 26) , 50 . 5° ( 11 . 52 ) , 55.7°(9.69) 

59.1° (7. 35). The asymptotic random intervals vary from -0.52 and 
+0.60 dB to -0.62 and +0.73 dB in Table 2. 


If the successive values of t R are independent in the 
probability sense, then the full table extension of Table 2 as 
the values vary slowly along the orbit could be used for all of 
the preceding analysis, and in particular Eqns. (4) , (22) , (24) , 
(27), (31), (34), (50), (51) and (52) would be theoretically cor- 
rect. 


8.4 Time correlated random variables . For a particular se- 
quence of 25 by 25 km areas a particular distance cross track 
along the swath, suppose that successive values of t R at the 
times tj, tj + At and so on as in, say, (60) are given by 

W’ t R (t l + At) » t R (t l + 2At) ’ t R (t l + nAt) 

It is possible for the successive values of t R to have a unit 

variance, zero mean normal PDF for a long enough sample, but the 
successive values can be highly correlated such that 

e(t R (t l ' 1 * t R (t l + nAt)) * ° 

and instead 

£ (t R (t 1 ) , t R (t 1 + nAt)) = p(nAt) 

where P is the correlation coefficient. 

The next simplest model of the error structure for R would 
then be a stationary Gaussian random process similar to the 
linear model for ocean waves sampled as a function of time at a 


(64) 


(65) 
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fixed point, but with a longer (?) , or different, time scale. 
For the four beams involved, the model might have to be extended 
to a vector Gaussian process. Very little is known about how 
this could be done, but the off nadir altimeter effects and the 
altimeter attitude data might provide guidance. 

It might be possible to investigate whether or not the er- 
rors in R are time correlated along the orbit once the scat- 
terometer is in orbit if the mesoscale variability in the wind 
vectors for each area sampled is not too great. The small scale 
variability of the wind will be described next. The analysis 
will be completed after this complication is described. 

8.5 Model errors, synoptic scales and mesoscales . The con- 
cept of a model for o° or as an extended concept the model for 
backscatter involving 0, X and 0 in Eqns. (7) and (8), is needed 
to recover vector winds from backscatter data. The area inside 
the 25 by 25 km area defining the array of the sampled 
backscatter values that is actually illuminated by the radar to 
obtain the backscatter estimate is much smaller than 625 km 2 , and 
the wind over this smaller area has a considerable amount of me- 
soscale variability from one 25 by 25 km area to the next. The 
concept of a synoptic scale perturbed from one area to the next 
by mesoscale wind fluctuations is useful. Although the larger 
ocean gravity waves may not change appreciably, the Bragg scat- 
terers will respond to these fluctuations. These concepts have 
been explored more fully by Pierson (1983a) , Pierson, et al. 
(1984) and summarized in Pierson, et al. (1986). Basically on an 
i, j grid of 25 by 25 km areas, the departures from a synoptic 
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scale wind Ug, X g could be represented by the pattern shown below 
near an arbitrary origin. 


5 S * AU 0,2’ X S* iX 0,2 

V AU 1,2- 

X S +AX 1,2 

°S +Ax 2,2 ; 

X S +AX 2 ,2 

V AU 0,1 ; X S +AX 0 , 1 

V Au i,r 

V AX 1,1 

°S +AU 2,1 J 

V AX 2,1 

°S * AU 0,0’ V AX 0,0 

V AU 1,0 ; 

x s +Ax i,o 

V AU 2,0 

; X S +AX 2 , 


The concurrent variability in the actual wind, the random 
variability of t p and the two possibilities (at least) for the 
variability of t D make it difficult to determine whether or not 
the different values of o° are independent random variables or 
correlated as a function of time and adjacent sampled areas of 
the ocean. 

Errors in the model are, at present, mostly systematic and 
not random. Whatever these errors are, they must be modeled in 
terms of o°, not 0° . Simply increasing VAR( o° ) by changing a , 6 
and y does not help because this changes the theoretical behavior 
of the radar and not the model. Errors in the theoretical model 
for backscatter from waves need to be treated some other way 
after the gross errors in the model for o° are removed. 
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8.6 A test for t constant over short orbit segments . 
Statistics can be computed from sample values of o° that may 
reveal whether or not successive values of o° are highly cor- 
related for some subset of the sample values. For a given an- 
tenna and row of 25 by 25 km areas parallel to the subsatellite 
track, the wind vectors may be slowly varying enough and the me- 
soscale variability may be small enough to verify the assumption 
that the values of t R are highly correlated. The failure of the 
test does not show that the t R are uncorrelated since mesoscale 
effects may mask the hypothesis. The test would then show that 
other sources of variability dominate even if the t R are cor- 
related and that the correlation can be neglected. 

For a sequence of n areas parallel to the subsatellite 
track, with n perhaps as large as 5 or 10, and for the same 
polarization and antenna, the available random variables would be 


a 


o 

V 




a 


o 

n 


Eqn. 


Suppose that t is constant 

K 

O 

(60) , the expected value of o 
e (a°) = + t R SD(o R ) 


for the n samples. Then from 
is given by Eqn. (66) . 


( 66 ) 


It follows that if 



then 


(67) 


^ < 5 ? = °M + t R SD( °R ) 


( 68 ) 
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and that 


e (o 


so that 


-o,2 

a P 


-2-^- VAR(o°) 
n p' 


(69) 


S 5 = rrrl^i-^' 


(70) 


has an expected value of VAR( ) and is an unbiased estimate of 

o 

VAR (Op ) . This ought to be close to 


= 0 ^( 5 °) 2 + o" ♦ y x 


(71) 


If, on the other hand, t and t are independent for this 

P K 


sample , then 


n 

-0 1 r O 

o ~ = — ) a . 

2 n J i 

is the sample mean such that 


(72) 


e ( 5 °) - 


(73) 


Also 


n 


C 2 1 V fn° «<>•/ 

S 2 ‘ ST^T | ( °i - °2> 


(74) 


is the unbiased sample variance. The quantity S 2 has an expected 
value equal to the variance of o° and is an unbiased estimate of 
VAR(c?°) which ought to be close to 


s 2 * = a 2 ^ 2 ^ + °2 + Y 2 


(75) 
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from Eqn. (62), 


The different samples of size n do not contain any informa- 
tion on whether either a ° , and or , s 2 and were sam- 
pled, but the statistics could contain this information. Only 
one sample mean is known and only one sample variance is known. 
The sample mean can be substituted into both Eqns. (67) and (71) . 

2 2 2 
If S 2 and S 2 * are close together, with perhaps S 2 larger 

than S 2 * then successive values of o° are probably independent, 

2 2 

and Table 2 is appropriate. If Sj and S ^ are close together and 

2 2 

if S 2 * is seven to twenty times larger than S then successive 

values of t R are highly correlated. Locally, the sampling 
variability is described by Table 1, but the expected values of 
the estimates do not equal the values for the model and can be 
biased either high or low. 

Thus, if the above statistics yield values for estimates of 
the variance perhaps 50% higher (or perhaps two to three times lar- 
ger) , but not seven to twenty times larger, than those in Table 1, 
the independence of successive values of t R is questionable. The 
use of Table 2 to recover the winds poses problems because the er- 
rors in the recovered winds from the MLE are systematic locally and 
not independent. This contingency needs to be treated theore- 
tically and modeled. A test of the hypothesis can be made rather 
quickly and the contingent plans implemented if needed. 

It is also possible that improved methods for reducing the ad- 
ditional sources of variability included in Table 2 may be 
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developed so that the values of a in that table may be reduced con- 
siderably prior to the launch of the spacecraft. All of the results 
presented prior to Part 8 are subject to revision if successive 
values of t R are not independent in the probability sense. 

8.7 The Amazon rain forest . The SEASAT-SASS antenna pattern 
was calibrated against homogenious and isotropic backscatter from 
the Amazon rain forest. These values would not vary because of me- 
soscale fluctuations in the wind over the ocean as described in Sec- 
tion 8.5. The test described in Section 8.6 can easily be applied 
over fairly long orbit segments so as to be able to determine whe- 
ther or not the added effects in Table 2, compared with Table 1, 
contribute to the variability of the sample values of the backscat- 


ter. 


I 

i 

I 

j 9. EXAMPLES BASED ON TABLE 2 FOR RANDOM INTERVALS AND 

RESULTS FOR A CANDIDATE MODEL FUNCTION 

; 9.1 Assumptions . In this part, it will be assumed, subject 

to revision, that Table 2 represents the parameters, a, g and y, 
for the various tabulated incidence angles. What could happen if 
Table 1 should have been used along with slowly varying values of 
t R will be mentioned at times. 

9.2 Random intervals . Tables 3.1 to 3.10 provide values 

for the end points of the random intervals with various 

o 

probabilities of including the model value, o M , given a sample 
o 

value of o • The calculations are based on the ten sets of 
values of a, 6 and y for the ten different incidence angles in 
Table 2. 

The first column is the value of t appropriate to the sample 

value, o°. The integral from t to 00 is the probability that 
o 

° M is zero (it cannot be negative) . 

The sample value of a 0 is in the second column. The column 
headed by 5% is the lower bound on a range of values for o° such 
that the value in that column and the value in the 95% column 
provide an interval that has a 90% chance of including the value 
of o°. 

The 15.87% and 84.13% columns are for l t 1 l = |t 2 l = !• The in- 
tervals defined by these end points have a 68.26% chance of en- 
closing the value of . 
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From Table 3.1 for the first row, a few of the statements 
that can be made are: 

(1) If o° = 0.18, the probability that = 0 is 

negligible. 

(2) P(0. 1155 < 0 ° < 0.2713) = 0.90 

(3) P ( 0 . 1387 < 0° < 0.2389) = 0.6826 

In contrast, if the sample value is 0.06 (i.e. o° = 0.06), 

o 

the probability that o M could be identically zero can be found in 
the tables for the normal PDF. Also 

(4) P(0.0077 < 0° < 0.0576) = 0.90 

M 

and so on. 

. . o 

If the sample value is 0.012, the probability that o M =0 is 
not negligible. Also 

(5) P(0 < a° < 0.0705) < 0.90 

and 

(6) P(0. 0705 < 0° < co) = o.05 

If 0 ° = -0.06, none of the intervals exist. The probability 
that = 0 is almost but not quite 1.00 as calculated from t and 
the probability that o° is greater than zero is less than 0.05. 

In terms of the analysis given in Section 5.4, the bottom 
four rows show the values of 0 ° such that, in order, the lower 
bounds on the 90% and 68.26% intervals become zero and the upper 
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TABLE 3.1 STATISTICS FOR RANDOM INTERVALS AND MAXIMUM 
LIKELIHOOD ESTIMATES (TIMES 10 2 ) 0 = 15.6 , 

a « 0.019125 B = 0.002517, Y - 0.0009917 


t J 


5% 

15.87% 

84.13% 

95% 

MLE 

5.72 

18.0 

11.55 

13.87 

23.09 

27.13 

17.54 

5.33 

16.8 

10.50 

12.77 

21.74 

25.65 

16.36 

4.95 

15.6 

9.44 

11.67 

20.40 

24.18 

15.19 

4.57 

14.4 

8.38 

10.56 

19.05 

22.71 

14.01 

4.19 

13.2 

7.32 

9.46 

17.71 

21.24 

12.83 

3.81 

12.0 

6.25 

8.35 

16.38 

19.78 

11.65 

3.43 

10.8 

5.17 

7.23 

15.05 

18.33 

10.48 

3.05 

9.6 

4.08 

6.11 

13.72 

16.87 

9.30 

2.67 

8.4 

2.98 

4.99 

12.40 

15.45 

8.12 

2.29 

7.2 

1.88 

3.86 

11.08 

14.02 

6.99 

1.91 

6.0 

0.77 

2.72 

9.77 

12.61 

5.76 

1.52 

4.8 


1.58 

8.46 

11.20 

4.59 

1.14 

3.6 


0.43 

7.16 

9.80 

3.41 

0.76 

2.4 

- 

- 

5.87 

8.41 

2.23 

0.38 

1.2 

- 

- 

4.59 

7.05 

1.06 

0 

0 


- 

3.31 

5.69 

0 

-0.38 

-1.2 

- 

- 

2.04 

4.35 

0 

-0.76 

-2.4 

- 

- 

0.78 

3.01 

0 

-1.14 

-3.6 


- 

- 

1.05 

0 

-1.53 

-4.8 

- 

- 

- 

0.41 

0 

-1.91 

-6.0 

— 

— 

— 

— 

0 

1.654 

5.18 

0 

NA 

NA 

NA 


1.000 

3.14 

NA 

0 

NA 

NA 


-1.000 

-3.14 

NA 

NA 

0 

NA 


-1.654 

-5.18 

NA 

NA 

NA 

0 



NA is "not applicable" 
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TABLE 3.2 STATISTICS FOR RANDOM INTERVALS AND MAXIMUM 
LIKELIHOOD ESTIMATES (TIMES 10 2 ) 6 = 19.9° , 

a = 0.01648, 8 = 0.0013575, Y = 0.0008225 


t 


5% 

15.87% 

84.13% 

95 % 

MLE 

6.29 

18.0 

12.22 

14.32 

22.42 

25.84 

17.64 

5.86 

16.8 

11.15 

13.21 

21.09 

24.40 

16.46 

5.44 

15.6 

10.08 

12.10 

19.76 

22.96 

15.28 

5.02 

14.4 

9.00 

10.98 

18.44 

21.53 

14.10 

4.60 

13.2 

7.91 

9.86 

17.12 

20.10 

12.92 

4.18 

12.0 

6.82 

8.73 

15.81 

18.67 

11.74 

3.77 

10.8 

5.72 

7.60 

14.50 

17.27 

10.56 

3.55 

9.6 

4.61 

6.47 

13.19 

15.87 

9.38 

2.93 

8.4 

3.49 

5.33 

11.89 

14.48 

8.20 

2.51 

7.2 

2.37 

4.19 

10.59 

13.09 

7.02 

2.09 

6.0 

1.23 

3.04 

9.30 

11.72 

5.84 

1.67 

4.8 

0.008 

1.88 

8.02 

10.35 

4.66 

1.26 

3.6 

- 

0.71 

6.74 

9.00 

3.48 

0.84 

2.4 

- 

- 

5.48 

7.66 

2.29 

0.42 

1.2 

- 

- 

4.21 

6.34 

1.12 

0.00 

0 

- 

- 

2.96 

5.02 

0 

-0.42 

-1.2 

- 

- 

1.72 

3.72 

0 

-0.84 

-2.4 

- 

- 

0.48 

2.44 

0 

-1.26 

-3.6 

- 

- 

- 

1.16 

0 

-1.67 

-4.8 

— 

— 

— 

— 

0 

1.654 

4.72 

0 

NA 

NA 

NA 


1.000 

2.87 

NA 

0 

NA 

NA 


-1.000 

-2.87 

NA 

NA 

0 

NA 


-1.654 

-4.72 

NA 

NA 

NA 

0 



NA is "not applicable" 
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TABLE 3.3 STATISTICS FOR RANDOM INTERVALS AND MAXIMUM 
LIKELIHOOD ESTIMATES (TIMES 102 ) 6 = 24.7 °, 

a = 0.01729, 8 = 0.0004497, Y = 0.00008354 


t 


5% 

15.87% 

84.13% 

95f 

MLE 

5.91 

5.40 

3.58 

4.24 

6.79 

7.87 

5.29 

5.51 

5.04 

3.26 

3.91 

6.39 

7.44 

4.93 

5.12 

4.68 

2.94 

3.58 

5.99 

7.01 

4.58 

4.72 

4.32 

2.62 

3.24 

5.60 

6.59 

4.22 

4.33 

3.96 

2.29 

2.91 

5.20 

6.15 

3.87 

3.94 

3.60 

1.96 

2.57 

4.80 

5.72 

3.52 

3.54 

3.24 

1.63 

2.23 

4.41 

5.29 

3.16 

3.15 

2.88 

1.30 

1.89 

4.01 

4.87 

2.81 

2.76 

2.52 

0.96 

1.55 

3.63 

4.45 

2.46 

2.36 

2.16 

0.63 

1.20 

3.23 

4.03 

2.10 

1.97 

1.80 

0.28 

0.86 

2.85 

3.62 

1.75 

1.58 

1.44 

- 

0.51 

2.47 

3.21 

1.39 

1.18 

1.08 

- 

0.16 

2.08 

2.80 

1.04 

0.79 

0.72 

- 

- 

1.70 

2.40 

0.69 

0.39 

0.36 

- 

- 

1.13 

1.80 

0.33 

0.00 

0 

- 

- 

0.95 

1.61 

0 

-0.39 

-0.36 

- 

- 

0.57 

1.21 

0 

-0.79 

-0.72 


- 

0.19 

0.83 

0 

-1.18 

-1.08 

- 

- 

- 

0.25 

0 

-1.58 

-1.44 

- 

- 

- 

0.07 

0 

-1.97 

-1.80 

— 

_ 

_ 

_ 

0 

1.654 

1.50 

0 

NA 

NA 

NA 


1.000 

0.91 

NA 

0 

NA 

NA 


-1.000 

-0.91 

NA 

NA 

0 

NA 


-1.654 

-1.50 

NA 

NA 

NA 

0 



NA is "not applicable" 
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TABLE 3.4 STATISTICS FOR RANDOM INTERVALS AND MAXIMUM 
LIKELIHOOD ESTIMATES (TIMES 10 3 ) 6 =31.4°, 

a = 0.01713 6 = 0.0001642, y = 9.405 *10-6 


t 


5% 

15.87% 

84.13% 

95% 

MLE 

5.87 

18.0 

11.89 

14.11 

22.68 

26.32 

17.62 

5.48 

16.8 

10.83 

13.00 

21.35 

24.87 

16.44 

5.09 

15.6 

9.75 

11.89 

20.02 

23.43 

15.26 

4.70 

14.4 

8.67 

10.77 

18.70 

21.99 

14.08 

4.30 

13.2 

7.59 

9.65 

17.37 

20.56 

12.90 

3.91 

12.0 

6.50 

8.52 

16.06 

19.14 

11.72 

3.52 

10.8 

5.40 

7.40 

14.74 

17.71 

10.54 

3.13 

9.6 

4.29 

6.27 

13.43 

16.31 

9.36 

2.74 

8.4 

3.17 

5.13 

12.13 

14.91 

8.18 

2.35 

7.2 

2.05 

3.99 

10.83 

13.52 

7.00 

1.96 

6.0 

0.91 

2.83 

9.54 

12.13 

5.87 

1.57 

4.8 

- 

1.60 

8.25 

10.77 

4 . 64 

1.17 

3.6 


0.52 

6.97 

9.41 

3.46 

0.78 

2.4 


- 

5.70 

8.06 

2.28 

0.39 

1.2 

- 

- 

4.44 

6.73 

1.10 

0.00 

0 


- 

3.18 

5.40 

0 

-0.39 

- 1.2 

- 

- 

1.93 

4.10 

0 

-0.78 

- 2.4 


- 

0.69 

2.80 

0 

-1.17 

- 3.6 

- 

- 

- 

1.52 

0 

-1.57 

- 4.8 

— 

— 

— 

0.26 

0 

1.654 

5.04 

0 

NA 

NA 

NA 


1.000 

3.07 

NA 

0 

NA 

NA 


-1.000 

-3.07 

NA 

NA 

0 

NA 


-1.654 

-5.04 

NA 

NA 

NA 

0 



NA is "not applicable". 
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TABLE 3.5 STATISTICS FOR RANDOM INTERVALS AND MAXIMUM 
LIKELIHOOD ESTIMATES (TIMES 10 3 ) 9 = 35.4°, 

a = 0.01775, 6 = 0.00011, Y = 2.478 * 10~ 6 


t 

SAMPLE a° 

5% 

15.87% 

84.13% 

95% 

MLE 

7.24 

12.0 

8.36 

9.66 

14.88 

17.17 

11.74 

6.76 

11.2 

7.66 

8.93 

13.99 

16.18 

10.95 

6.27 

10.4 

6.96 

8.20 

13.09 

15.20 

10.17 

5.79 

9.6 

6.26 

7.46 

12.19 

14.22 

9.38 

5.31 

8.8 

5.55 

6.73 

11.30 

13.25 

8.59 

4.83 

8.0 

4.84 

5.99 

10.41 

12.28 

7.81 

4.34 

7.2 

4.12 

5.25 

9.53 

11.32 

7.02 

3.86 

6.4 

3.40 

4.50 

8.64 

10.36 

6.24 

3.38 

5.6 

2.67 

3.75 

7.76 

9.41 

5.45 

2.90 

4.8 

1.94 

3.00 

6.88 

8.46 

4.66 

2.41 

4.0 

1.20 

2.24 

6.01 

7.52 

3.88 

1.93 

3.2 

0.45 

1.48 

5.16 

6.59 

3.09 

1.45 

2.4 

- 

0.72 

4.28 

5.67 

2.30 

0.97 

1.6 

- 

- 

3.43 

4.75 

1.52 

0.48 

0.8 

- 

- 

2.57 

3.85 

0.73 

0.00 

0 

- 

- 

1.73 

2.96 

0 

-0.48 

-0.8 

- 

- 

0.89 

2.07 

0 

-0.97 

-1.6 

- 

- 

0.06 

1.20 

0 

-1.45 

-2.4 

— 


— 

0.35 

0 

1.654 

2.73 

0 

NA 

NA 

NA 


1.000 

1.66 

NA 

0 

NA 

NA 


-1.000 

-1.66 

NA 

NA 

0 

NA 


-1.654 

-2.73 

NA 

NA 

NA 

0 



NA is "not applicable". 
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TABLE 3.6 STATISTICS FOR RANDOM INTERVALS AND MAXIMUM 
LIKELIHOOD ESTIMATES (TIMES J.0 4 ) 6 = 41°, 

a = 0.01841, B = 5.53 X 10 , Y = 5.849*10' 


t 


5% 

15.87% 

84.13% 

95% 

MLE 

6.28 

48.0 

31.97 

37.73 

60.64 

70.64 

46.87 

5.86 

44.8 

29.17 

34.80 

57.05 

66.70 

43.73 

5.44 

41.6 

26.36 

31.86 

53.46 

62.78 

40.58 

5.02 

38.4 

23.53 

28.92 

49.88 

58.87 

37.44 

4.60 

35.2 

20.69 

25.97 

46.32 

54.98 

34.30 

4.18 

32.0 

17.83 

23.00 

42.76 

51.10 

31.16 

3.77 

28.8 

14.95 

20.03 

39.22 

47.24 

28.01 

3.35 

25.6 

12.05 

17.04 

35.68 

43.41 

24.87 

2.93 

22.4 

9.13 

14.04 

32.16 

39.59 

21.73 

2.51 

19.2 

6.18 

11.03 

28.66 

35.80 

18.58 

2.09 

16.0 

3.21 

7.99 

25.17 

32.04 

15.44 

1.67 

12.8 

0.21 

4.95 

21.70 

28.31 

12.30 

1.26 

9.6 

- 

1.88 

18.24 

24.61 

9.16 

0.84 

6.4 

- 

- 

14.81 

20.94 

6.02 

0.42 

3.2 

- 

- 

11.40 

17.31 

2.88 

0.00 

0 

- 

- 

8.00 

13.72 

0 

-0.42 

-3.2 

- 

- 

4.64 

10.17 

0 

-0.84 

-6.4 

- 

- 

1.30 

6.66 

0 

-1.26 

-9.6 

— 

— 

— 

3.19 

0 

1.654 

12.58 

0 

NA 

NA 

NA 


1.000 

7.65 

NA 

0 

NA 

NA 


- 1.000 

-7.65 

NA 

NA 

0 

NA 


-1.654 

-12.58 

NA 

NA 

NA 

0 



NA is "not applicable". 
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TABLE 3.7 STATISTICS FOR RANDOM INTERVALS AND MAXIMUM 
LIKELIHOOD ESTIMATES (TIMES 10 3 ) 0 = 46° , 

a = 0.019375, 6 = 0.000043, T = 2.688 X 10 


t ; 


5% 

15.87% 

84.13% 

95% 

MLE 

5.79 

3.0 

1.93 

2.31 

3.85 

4.53 

2.92 

4.82 

2.5 

1.49 

1.86 

3.29 

3.91 

2.43 

3.86 

2.0 

1.05 

1.39 

2.73 

3.30 

1.94 

2.89 

1.5 

0.597 

0.929 

2.17 

2.69 

1.45 

1.93 

1.0 

0.137 

0.459 

1.62 

2.10 

0.960 

1.157 

0.6 

- 

0.078 

1.19 

1.63 

0.568 

0.965 

0.5 

- 

- 

1.08 

1.51 

0.4695 

0 

0 

- 

- 

0.546 

0.940 

0 

-0.965 

-0.5 

- 


0.002 

0.383 

0 

-1.544 

-0.8 

— 

— 

— 

0.057 

0 

1.645 

0.853 

0 

NA 

NA 

NA 


1.000 

0.518 

NA 

0 

NA 

NA 


-1.000 

-0.518 

NA 

NA 

0 

NA 


-1.645 

-0.853 

NA 

NA 

NA 

0 



NA is "not applicable". 
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TABLE 3.8 STATISTICS FOR RANDOM INTERVALS AND MAXIMUM 

LIKELIHOOD ESTIMATES (TIMES 10 3 ) 6 = 50.5°, 7 

a = 0.019917, 3 = 0.000031, Y = 1.539 X 10 


t 

SAMPLE o° 

5% 

15.87% 

84.13% 

95% 

MLE 

6.12 

2.4 

1.57 

1.87 

3.06 

3.59 

2.34 

5.10 

2.0 

1.23 

1.50 

2.61 

3.09 

1.95 

4.08 

1.6 

0.872 

1.14 

2.16 

2.60 

1.55 

3.06 

1.2 

0.512 

0.765 

1.72 

2.11 

1.16 

2.04 

0.8 

0.144 

0.389 

1.28 

1.64 

0.769 

1.02 

0.4 

- 

0.007 

0.840 

1.17 

0.377 

0 

0 

- 

- 

0.412 

0.709 

0 

-1.02 

-0.4 

— 

— 

— 

0.265 

0 

1.645 

0.645 

0 

NA 

NA 

NA 


1.000 

0.392 

NA 

0 

NA 

NA 


-1.000 

-0.392 

NA 

NA 

0 

NA 


-1.645 

-0.645 

NA 

NA 

NA 

0 



NA is "not applicable". 
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TABLE 3.9 STATISTICS FOR RANDOM INTERVALS AND MAXIMUM 

LIKELIHOOD ESTIMATES (TIMES 10 3 ) 6 =55.7°, _ 7 

a = 0.021640 @ = 0.000037595, Y - 1.709 X 10 



5% 

15.87% 

84.13% 

~95l 

MLE 

5.81 

2.4 

1.53 

1.84 

3.10 

3.68 

2.33 

4.84 

2.0 

1.18 

1.48 

2.65 

3.17 

1.94 

3.87 

1.6 

0.833 

1.11 

2.20 

2.67 

1.55 

2.90 

1.2 

0.476 

0.741 

1.75 

2.18 

1.16 

1.94 

0.8 

0.111 

0.367 

1.31 

1.70 

0.765 

1.161 

0.48 

- 

0.064 

0.956 

1.31 

0.452 

0.967 

0.4 

- 

- 

0.869 

1.22 

0.373 

0 

0 


- 

0.438 

0.757 

0 

-0.967 

-0.4 

- 

- 

0.0140 

0.307 

0 

-1.549 

-0.64 

— 

— 

— 

0.0434 

0 

1.645 

0.680 

0 

NA 

NA 

NA 


1.000 

0.413 

NA 

0 

NA 

NA 


-1.000 

-0.413 

NA 

NA 

0 

NA 


-1.645 

-0.680 

NA 

NA 

NA 

O 



NA is "not applicable". 
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TABLE 3.10 STATISTICS FOR RANDOM INTERVALS AND MAXIMUM 

LIKELIHOOD ESTIMATES (TIMES 10^) 6 = 59.1°, _ 7 

a = 0.02376, 6 = 0.00005929, Y = 3.055 X 10 


t 


5% 

15.87% 

84.13% 

95% 

MLE 

5.43 

3 

1.84 

2.25 

3.95 

4.74 

2.90 

4.52 

2.5 

1.41 

1.80 

3.38 

4.10 

2.41 

3.62 

2.0 

0.98 

1.35 

2.81 

3.47 

1.93 

2.71 

1.5 

0.535 

0.886 

2.25 

2.84 

1.44 

1.81 

1.0 

0.083 

0.421 

1.69 

2.23 

0.948 

1.085 

0.6 

- 

0.045 

1.25 

1.74 

0.557 

0.905 

0.5 


- 

1.14 

1.62 

0.460 

0 

0 

- 

- 

0.59 

1.02 

0 

-0.905 

-0.5 

- 

- 

0.056 

0.456 

0 

-1.629 

-0.9 

— 


— 

0.010 

0 

1.645 

0.909 

0 

NA 

NA 

NA 


1.000 

0.553 

NA 

0 

NA 

NA 


-1.000 

-0.553 

NA 

NA 

0 

NA 


-1.645 

-0.909 

NA 

NA 

NA 

0 



NA is "not applicable". 
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bounds on the 68.26% and 90% intervals become zero. All of the 
probability statements in Section 5.4 can be made for the various 
values of the random variable, o°, based on statistics calculated 
from this random variable as derived in Section 7. 

The last column in each table provides the values of the max- 
imum likelihood estimates for each sample value of o°. Each 

estimate is slightly less than the value of o° . if the estimate, 

o . . . o 

o , is negative, the maximum likelihood estimate of a M must be 

zero. 

There is a theoretical inconsistency in the definition of 
the maximum likelihood estimate as in Eqn. (24) . If there is a 
finite probability that a° = 0, since the PDF is otherwise con- 
tinuous, the MLE strictly defined ought to be zero. In another 
sense the PDF is maximized for two different value of o° and the 
realistic choice is the one tabulated. 

Figs. 4a to 4g are graphs of these values for incidence 
angles of 15.6°, 19.9°, 24.7°, 31.4°, 41.0°, 50.5° and 59.1°. 

The graphs are entered with the sample value of a° on the ver- 
tical scale and the curves then define the various random inter- 
vals with the indicated probabilities of enclosing the value of 
as given on the horizontal axis. Supplementary values in dB 
are also shown for each graph. 

For Fig. 4a, the sample values range from 0.20 (-6.99 dB to) 

zero (- 00 dB) to -0.10 (dB not defined). If c° has the value 

% 

given by 1.654 Y , the 90% random interval has a range from - 00 to 
somewhat greater than -10 dB. Trying to understand backscatter 
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SAMPLE 


VALUE 




68.267. 


907. 


EXPECTED VALUE 


0.10 0.15 0.20 0.25 


(-10.0) (-8.24) (-6.99) (-5.23) dB 


6 = 15.6 


0*0.019125, £*0.002517, y *0.0009917 


osymptotic volues t*7.23, K«*0.I38 


dB -0.56, +0.65 


SAMPLE 


VALUE 


•v/ 




68.267. 


EXPECTED VALUE 


0.20 0.25 


B * 19.9 


= 0.01648, £=0.0013575, y =0.0008225 
osymptotic volues t = 7.79 K~ = 0.128 
B -0.52, +0.60 


rrrrs 



































Fig. 4c Same as Fig. 4a Except for 6 = 24.7°. 



Fig. 4d Same as Fig. 4a Except for 6 = 31.4°. 


61 









Fig. 4e Sane as Fig. 4a Except for e * 41.0° 



Fig. 4f Sane as Fig. 4a Except for 6 *= 50.5° 
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Fig. 4g Same as Fig. 4a Except for 6 = 59.1°. 
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values in terms of decibels or logarithms of o° and o° can be dif- 
ficult at times. 

These seven graphs at first glance appear to be very much 
alike, but closer inspection shows them to be quite different. 
The ranges of a° and a° are different. The intercepts with the 

M 

0 • • \* • 
a axis, which depend on y 2 are different and the lengths of the 

random intervals increase at different rates depending on the 

size of a . 

Since y decreases by only 1.5 to 2% when Table 1 is compared 
to Table 2, had Table 1 been used to construct the graphs, the 
four curves would have had slopes closer to unity and the lengths 
of the random intervals would not increase as rapidly as they do 
in these figures. The variability introduced by the uncertainty 
in the quantity, R, is thus important for high values of a 0 (and 
o° ) , in a relative sense, as determined especially by a, and by 
3 and y . 


9.3 Interpretation in terms of a candidate model . For il- 
lustrative purposes, the candidate model will be the one given by 
Donelan and Pierson (1987). The curves for the model are from 
Fig. 7 of this reference. Figs. 5a and 5b are for a water tem- 
perature of 0°C and incidence angles respectively of 20° , 40° and 
60° and 30° and 50°. Figs. 5c and 5d are for a water temperature 
of 30° and the same incidence angles as in the previous two 
figures. The slight differences in incidence angle between Table 
2 and these values have been neglected. All four of these 
figures are based on the sampling variability of one single 25 by 
25 km area for vertical polarization. 
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Figs. 6a to 6d show the result of averaging the a° values 
from four 25 by 25 km areas in a row parallel to the subsatellite 
track. The four incidence angles and the four values of a , 3 and 
Y would be averaged also. The values of a, 3 and Y have been 
divided by four based on (27) . 

For purposes of discussion, it will be assumed that the 
model is correct. The results are for either one 25 by 25 km area 
at the indicated incidence angle and temperature as in Figs. 5a 
to 5d or for a row of four areas parallel to the subsatellite 
track as in Figs. 6a to 6d. The values of a, 3 and y from Table 2 
have been used for the incidence angles closest to those in these 
figures. 

For each incidence angle in Figs. 5a to 5d there are five 
curves. The center curve is for the model, and it is surrounded 
by the four curves labled 5%, 15.83%, 84.13% and 95%. The 
backscatter estimate is assumed to be for an upwind look. The con- 
clusions of Donelan and Pierson that the backscatter can become 
undetectable before the wind becomes zero, that there is no power 
law and that backscatter may saturate and then decrease for high 
enough winds are illustrated by this set of four figures. 

The interpretation of each one of the eight figures depends 
on the meaning of the vertical axis. If for upwind repre- 
sents the model value, then entering with the wind speed on the 
horizontal axis provides the model value for on the vertical 
axis. Eqn. (12) can then be used to calculate the scatter of the 
random variable o° about the value of and graph the four 
curves with the assigned probabilities for Figs. 5a to 5d. Thus 
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for Fig. 5a if the wind is, say, 10 m/s for an incidence angle of 

o o 

40 , the random variable, a will be between about -14 dB and -17 
dB 90% of the time. For Figs. 6a to 6d, Eqn. (27) yields a 
reduced variance and standard deviation, and this new value is 
used in (12) to calculate the 5% and 95% curves. At 10 m/s for, 
example, in Fig. 8a at a 40° incidence angle, the random variable 
will be within a range from about -15.4 dB to about -16.8 dB 90% 
of the time. 

For all eight figures for low enough wind speeds, there is 
no backscatter predicted by the model above -47.5 dB. In the 
model being used as an example, the prediction is that there is no 
backscatter from the waves for winds below a certain water tem- 
perature dependent wind speed. Below this threshold wind speed 
the model value for o° is zero. In Figs. 5a to 5d, the 5% and 
15.87% curves are drawn for 1 m/s at the values 1.645Y 2 and Y 2 

\s 

in dB. In Figs. 6a to 6d, the 5% curve is at 1.645 (y/ 4) 2 and is 
thus 3 dB lower for low winds than the corresponding curve in 
Figs. 5a to 5d. 

The 84.13% and 95% curves in the various Figs. 5a to 5d go 
to minus infinity dB at those points on the 0, axis in the ap- 
propriate Figure from Figs. 4a to 4g where these curves intersect 
the axis for one 25 by 25 km area. For low winds with a° in 
antilog form, the 90% range of values of o° would be from the ap- 
propriate antilog in the figure as an upper bound to its negative 
value as a lower bound. 

The lower left region for each incidence angle has the 
property that positive sample values of the backscatter from ac- 
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tual data plotted in dB against meteorologically obtained winds 
would scatter all over the region. This feature of actual data 
is well illustrated in Schroeder, et al. (1982). Negative values 
were discarded. 

For the range of wind speeds that were used for the study of 
the SEASAT SASS data, it would be possible to draw a number of 
different power laws within the 5% to 95% curves for Figs. 5a to 
5d. Graphs such as these figures about model lines as power laws 
would still form the same patterns for the low wind, low a 0 
areas, and thus light winds would still be uncertain, but well 
within present design requirements, except, perhap, for cold 
water. 

For one 25 by 25 km area (1 cell) , Figs. 5a to 5d show that 
if the incidence angle and water temperature are known, 90% of 
the sample values of a 0 will occur somewhere in the region below 
the 5% curve and above, or to the left, of the 95% curve. Were 
both the value of o° and the wind speed known the plotted point 
(0, o° ) would usually not fall on the model curve for a°. 

For four cells in a row parallel to the subsatellite track, 
the advantage of an average (if the wind does not change too 

much) is evident in Figs. 6a to 6d. Ninety percent of the points 

-o 

(U, o ) scatter in a much narrower range about the model curve 
and the area in the lower left for each incidence angle is made 
smaller. 

Figures similar to these eight figures could also have been 
prepared for averages of two c° values parallel to the subsatel- 
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lite track. The values of a, g and y would be reduced by a factor 
of two. The asymptotic values in Table 2 would all be reduced. 
The value of y 2 would also be reduced. As an example, for, say a 
35.4° incidence angle, instead of -0.54 and +0.62 dB, the asymp- 
totic values would be -0.39 and +0.43 and y 2 in dB would be -29.3 
instead of -27.8, dB. 

Each of these figures for a particular incidence angle could 
be extended to a three variable plot of o° + |t 2 |SD(a°), a° an( j 
°M _ I t l l SD *- a °p versus U and X. For a fixed U, the resulting sur- 
face for o° would vary as sketched in Figs. 5 and 6 and the upper 
and lower ranges for the scatter in o° would follow the undula- 
tions of the surface. 

If the vertical scale on these eight figures represents the 
sample value, either o° or o°, there is the temptation to enter 
one of these figures, with the sample value, find the values for 
the 5% and 95% curves, for example, and the corresponding wind 
speeds and conclude that, if looking upwind, there is a 90% 
chance that the wind speed is between these two values. As a 

first approximation, this proceedure yields a range of wind 
speeds that is not quite correct because the slopes of the curves 
are continuously changing with increasing wind speed. 

Figs. 7a to 7d illustrate a correct way to interpret these 
figures given a sample value of either o° or o° . Since the in- 
cidence angle is known, the values of a , 6 and y are known. The 
random interval with, say, a 90% chance of enclosing the value, 
o° can be computed, and the values for the end points will inter- 
sect the curve for o° in the various ways sketched in the 
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Fig. 7a) An Example of the Use of a 90% Random Interval for a 

Portion of one of the Sets of Curves in Either Fig. 5 

• or 6. Not to Scale; 

b) The Use of a 90% Random Interval if the Model Values 
Saturate and Then Decrease; 

c) The Use of a 90% Random Interval Near the Roll-Over 

Portion of a Model, the 5% and 95% Curves Have No 

Mscinin^ j 

d) A Possible 90% Random Interval Near the Maximum of the 
Roll-Over for a Model, the 5% and 95% Curves Have No 
Meaning. 
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figures, if the lower bound lies above the 5% curve for the 
chosen incidence angle and sample value. 

Fig. 7a illustrates the condition that the sample value is 
low enough so that the possible effect of saturation and roll- 
over for the model need not be considered. The upper and lower 
bounds on the random interval with a 90% chance of containing the 
value, a° , intersect the model curve at two points and the pos- 
sible values of a°, 90% of the time, will be on the heavy portion 
of the curve shown in Fig. 7a. Compared with the points where 
the line for o° intersects the 5% and 95% curves, the bounds on 
the range of wind speeds are broader. 

Should the effect of roll-over be even more pronounced than 
modeled in Figs. 5 and 6, three possibilities are shown in Figs. 
7b, 7c and 7d. For Fig. 7b, there are two possible ranges of 
values for a° and two possible ranges of wind speed. For Fig. 
7c, a very broad range of wind speeds is possible near the max- 
imum of the model curve. 

Fig. 7d is an interesting result. The sample value can be 
higher than the model value such that only the lower bound inter- 
sects the a° curve. If the model is correct, the range of wind 
speeds can be much smaller compared with Fig. 7c. 

The lower left corner for each incidence angle in Figs. 5 
and 6, is difficult to interpret for low winds and low 
backscatter estimates because of the distortion introduced by the 
logarithmic scales. Figs. 8a to 8d show how the sample values 
can be used for linear scales. For Figs. 8a to 8d, the main part 


74 


of the figure is the same. The curves for o° equal to o° and the 
5% and 95% bounds are sketched, but would need to be calculated 
carefully. Only the values of o° and the lines for the end points 
of the random interval with, say, a probability of 90% of en- 
closing the value, 0 °, change. 

For the model being used for an example, the value of 0 ° 
goes to zero at some water temperature dependent wind speed, de- 
signated by U T * For wind speeds less than U T , since there is no 
backscatter from the waves, the random variability is determined 
solely by the value of y in Table 2. Ninety percent of the sample 
values will scatter somewhere in the rectangular area bounded by 
0 < U < U T and + 1 . 645 y 2 and - 1.645y 2 if U is less than U T . 

For Fig. 8a, the sample value, o°, lies just above this rec- 
tangular area, but the lower bound is within this area. Never- 
theless, the values for the upper and lower bounds on the random 
interval intersect the curve for 0 ° at two points corresponding 
to speeds greater than U T . The range of wind speed for a 90% in- 
terval would be quite small. 

For Fig. 8b, the sample value is below the 5% line, the 
lower bound for o° is zero and the upper bound for a° exists. 
The 90% range on U is from zero to this upper bound. 

In Fig. 8c, the sample value is zero. An upper bound for o° 
exists. The range of wind speeds is from zero to the speed deter- 
mined by the upper bound. 

In Fig. 8d, the sample value is negative. The upper bound 
still exists and the wind speed range is from zero to this upper 
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Fig. 8a) The Use of the 90% Random Interval With Linear Scales 

for Values Near the Water Temperature Dependent Thresh- 
old Wind for a Sample Value Greater Than the 5% Curve; 

b) Same as a Except for a Positive Sample Value Below the 
5% Curve? 

c) Same as a Except for a Sample Value of Zero. The 
Points Marked A and B Correspond to the Intercepts 
of the Lower and Upper Bounds with the o° and o M 
Axes of Fig. 4; 

d) Same as a Except for a Negative Sample Value. 

For b, c and d, the Interval for o® is From Zero to 
the Upper Bound and the Probability 3 is 95%. 
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bound . 


9.4 The effect of constant orbit errors for fairly long 
times . If the statistical test given in Part 8.6 shows that the 
values of t are sensibly constant over a fairly large portion of 
a sampled area of the ocean. Figs. 5 and 6 would not be ap- 
plicable. There would be a curve on each of these figures for 
each incidence angle somewhat above or somewhat below the curve 
for o° (perhaps similar to a curve for o° for a slightly dif- 
ferent incidence angle) that would correspond to Eqn. (75) . The 
sample values would scatter about this curve with standard devia- 
tions calculated from Table 1 for the area of the ocean being sam- 
pled. The values of t R need not be the same for the eight dif- 
ferent antennas on the spacecraft. The analysis carried out 
above could perhaps be modified to account for this effect. The 
overall outcome for this possibility would be that wind 
recoveries could be systematically biased from one area of the 
ocean to another, especially for moderate and high winds. 
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10. MAXIMUM LIKELIHOOD ESTIMATES FOR 0 AND X . 


o 

10.1 Wind speed and direction conventions . Let the 45 an- 
tenna in Fig. 2 for vertical polarization be pointing at an angle 

of X clockwise relative to north. The values of X in (7) and 

o 

(8) are to be added to X Q to obtain the wind direction relative 
to north. The result is the direction from which the wind is 
blowing in degrees measured clockwise from north. For the two 
beams at 45° , in Fig. 2 the wind direction is simply X Q + X if X 
is known. 

When the beam at 115° to the subsatellite track samples the 
25 by 25 km area a few minutes later, if the wind direction rela- 
tive to the 45° beam was 0°, then the wind direction relative to 

the 115° beam would be -70° or 290°. Similarly, for the 135° 

o 

beam, the direction relative to this beam would be -90 . It is 
important that these sign conventions be correctly accounted for 
in the interpretation of the data. A transformation of variables 
to refer all directions to the 45° beam is accomplished by using 
X -70° and X -90° as the case may be in (7) or (8). 

10.2 Theories and models . There have been numerous efforts 
to apply electromagnetic theory combined with various properties 
of the waves on the ocean to the prediction of the backscatter 
that would be estimated by a scatterometer , or any other properly 
designed and calibrated radar. The problem is very complex since 
there are many physical processes involved, some of which may not 
even have been considered. 
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So as to avoid still more complicated notation, subscripts 
will be added to (7) and (8) to distinguish among the four dif- 
ferent measurements. A theory, or a model, essentially attempts 
to predict the quantity o° as a function of 0 , X and 6 as defined 
in Section 2.1. 


To determine the wind speed and direction, the information 
available for a 25 by 25 km area can be summarized as follows: 

a yyi is the estimate for vertical polarization from the 45° 
forward beam with ©j , , 3j and Yj known and a predicted ex- 
pected value from the theory, or model, given by X * V 

which should equal the expected value of this estimate. 

o 

a HH2 the estimate for horizontal polarization from the 

45° forward beam with © 2 , a 2 , 3 2 and y 2 known and a predicted 

expected value from the theory, or model, given by °hh'°- x, e 2 ) 
which should equal the expected value of this estimate. 


a HH3 is the estimate for horizontal polarization from the 
115 aft beam with 3^ an< * Y 3 known and a predicted ex- 


pected value from the theory, or model, given by 
°HH^ > X “ » ® 3 ) which should equal the expected value of this 
estimate. 

a° V4 is the estimate for vertical polarization from the 135° 
aft beam with 0 4 , a 4 , 3 4 and Y 4 known and a predicted expected 
value from the theory, or model, given by o° vv (0, X-90°, 0 4 ) which 
should equal the expected value of this estimate. 
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Each of the four expected values for the model backscatter 
is different because for the same wind speed, the model varies as 
a function of 6 and X, for the same X, the model varies as a func- 
tion of wind speed and 0 and for the same 0 the model varies as a 
function of wind speed and X. (There may be other physical quan- 
tities involved whose effects have already been accounted for.) 
Since however 6 is known, a known value of X will determine a 
unique wind speed from the model, and a known value of U will de- 
termine four, or fewer, unique values of the wind direction from 
the model even for just one of the estimates. 

10.3 Maximum likelihood estimates of 0 and X . Since the 
expected values and variances of the four different estimate of 
the backscatter are all different, the four different estimates, 

°VV1 ' °HH2 ' °HH3 and °VV4 are each a sain P le of size one from 
populations defined by (4) such that the expected values and 

variances are all different. 

The joint probability density function for the four sample 
values is the product of four probability density functions 
similar to (4) but with different expected values and variances. 
These depend on different known values for 0, a, 3 and Y as in 
Eqn. (76) where the expected values and variances are now 
predicted by the theory, or model. The four sample values are 
independent in the probability sense. 



o o 
°HH2’ a HH2 ' 


0 

a HH3’ 


f(Oyyj)* f(CT uuo ) ‘fXCTjju,) *f(0 Wyt ) 


VV4 


HH2 ; 


HH3 


W4 J 


(76) 


80 



Since the expected values for the four backscatter estimates 
are now given by a model, the maximum likelihood estimates of ii 
and X can be found by maximizing the right hand side of Eqn. (76) 
by varying U and X because o° vl , o° H2 , oj^ and o° ry4 are the 
known random variables. The steps are similar to those for Eqn. 
(24) except that the model, or theory, has imposed a further re- 
striction on the expected values and variances, or standard 
deviations, that can be used to maximize the probability. 

The result is given by Eqn. (77) as U and X are varied. 


MLE(U, X) = MAX - 


r 


C°Wl' X* 6j)) 

VARCoJyfO, X, 6j)) 


(0 HH2 ~ °HH (0> X> 6 2^ 2 

VARfo^fO, X, e 2 )) 


+ ^°HH3 ~ °HH^* X ~ 70 * e 3^ + 

VAR(o2 H (U, X-70°, e 3 )) 


( °VV4 - °W (0 > X - 90 °» 6 4 ))2 
VARCoJytO, X-90°, e 4 )) 


4 in ((VARCoJyfO, X, 6 1 ))/(a 1 (oJ vl ) 2 ♦ Bj 0^ ♦ Yj))) 

+ *» ((VARtoJ^CU, X, e 2 ))/(ct 2 (0 HH2 )2 + 6 2 °HH2 + Y 2 ))5 
♦ *> CtVABto^CU, x-70°. e 3 )/(a 3 (4 s ) 2 . 6 3 oS H3 . y 3 ») 
♦*» ((V**(o^,(0, X-90°, e 4 )/(a 4 (aj v4 ) 2 ♦ 6 o^, 4 . y 4 ») 


( 77 ) 
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The one half in Eqn. (24) has been omitted above because it is 
irrelevant. It could also have been omitted in Eqn. (24) . 


Eqn. (77) can have more than one local maximum. A computer 
search for these maxima can find them quickly. The power law as- 
sumption that was used for the model for backscatter during the 

SEASAT-SASS program and the use of only two backscatter values 

o o 

most of the time, such as and °yy 4 with the other terms 

omitted in (77), yields four maxima for many sample values corres- 
ponding to four different wind vectors such as Uj , Xj ; U 2 , X 2 ; 

0_; X_ and 0,, X, with wind vectors, one per quadrant, rela- 

tive to X q . The subject is discussed in considerable detail in 
Pierson (1983b) . The design in Figs. 1 and 2 is a simpler one 
than the one analysed by Pierson and Salfi (1982) when the NOSS 
(National Oceanographic Satellite System) was under con- 
sideration. It is not essential to know the inverse of, say. 




(0, X, 6) 


(78) 


in terms of 

0 = 0(0^, X, 0) (79) 

analytically in order to determine the maximum likelihood 
estimates of 0 and X. 

The maximum of the local maxima for Eqn. (77) may provide a 
unique wind vector most of the time. The beam at Ilf? for the 
example chosen above will almost always eliminate two of the four 
"X" type SASS solution. Differences betweeen horizontally 
polarized backscatter and vertically polarized backscatter could 
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eliminate one of the remaining two. 


It is also possible to extend these results to, say, four 25 
by 25 km areas that form a 50 by 50 km area. Eqn. (77) would 
then contain four times as many terms, one for each 25 by 25 km 
area. Averaging, as described above, could help reduce the 
putational load. 


com- 


11 SYSTEM AND MODEL VALIDATION 


11.1 Comparison with conventional wind measurements . Once 
the MLE has obtained values of (j and x» the resulting wind vector 
can be compared with a wind measured by conventional means in 
terms of an effective neutral wind computed from the conventional 
measurement for an assigned height, say 19.5 or 20 m, and the 0 
from a boundary layer model for the same height. The usual 
procedure is to make scatter plots of the two different speeds 
and the two different directions. It was difficult to obtain com- 
parison data for high winds for the SEASAT-SASS program. Also 
certain kinds of systematic errors are not immediately obvious in 
such analyses. 

The system and model validation for some future scat- 
terometer may be easier than it was for the SEASAT-SASS because 
of the expanded network of deep ocean data buoys deployed by the 
National Data Buoy Center. It is planned to have these data 
buoys obtain five consecutive 10 min wind averages followed by an 
8 min average to allow for data transmission each hour. Scat- 
terometer data can be co-located with a data buoy in both time 
and space, and the data buoy winds can be averaged for as long a 
time as appropriate to obtain a more representative wind, if 
needed. Gilhousen (1986) has reported on the present accuracy of 
wind measurements by the National Data Buoys. 
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There are plans to try to modify 18 National Data Buoys 
during the next few years so as to accomplish this change (Per- 
sonal communication, J. Wilkerson 1987) . With two 600 km wide 
swaths perhaps as many as 80% of the buoys will be in the swath 
twice each day so that in about one year there will be over 9000 
data buoy reports to be compared with SCAT recoveries. The 
various incidence angles will range over about 25° so that 360 
comparisons per incidence angle should be possible. 

Nevertheless, scatter plots of O(SCAT) versus U(BUOY) and 
X(SCAT) versus X(BUOY) will produce points that do not lie on the 
line of perfect fit. The reasons for this are many such as: 


(1) The model, or theory, is wrong. 

(2) The values of o° are estimates and do not equal the 
expected value predicted by a correct model. 

(3) The conventional measurements have systematic and 
sampling variability "errors". 


For the final validation of any model, it is essential to be 
able to compare the winds recovered from the backscatter 
estimates with the winds measured by a conventional system. The 
sources of "bias" and errors must be identified, and the effects 
of sampling variability for both the scatterometer winds and the 
conventional winds must be accounted for separately. 


The maximum likelihood estimates of 0 and X when compared 
with the 0 and x obtained from a data buoy provide not only these 
comparisons but also ways to compare the backscatter estimates 
from the scatterometer with the model backscatter values 
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recovered from the maximum likelihood estimate and the 

backscatter values that would be computed from the model from the 
data buoy winds with due consideration of the "errors" and sam- 
pling variability of the data buoy winds. 

For a given 25 by 25 km area that contains a data buoy, and 

for other areas that surround it, the available data (with the 
incidence angle known) consist of the following: 

Og, x B from the buoy. 

0 , X from the maximum likelihood estimate. 

°W'^B ’ ’ °HH^ ^B’ ’ °w^B’ computed from 

the model and the buoy wind for the appropriate wind 
direction relative to each beam. 

0^.(0, X), o° H (0, X), 0^(0, x). 0° V (U, X) from the wind speed 

and direction recovered from the maximum likelihood 
estimate. 

°\rvi’ °HH2’ a HH3’ °VV4 the randoin variables obtained by the 
scatterometer for the area. 

To avoid complicated subscripts, the values are ordered in 
the same sense as the first four terms in Eqn. (77) . The data 
for a large number of data buoy reports can be stratified ac- 
cording to and 6 for comparison purposes as the data set in- 

creases. Other buoy data such as wave spectral data, air tempera- 
ture and water temperature will also be needed. 

In general, Ug, X g will not agree exactly with 0, X from the 
maximum likelihood estimate, nor will the triplets of backscatter 
values for a given scatterometer beam such as ’ 

°vv(0, and °vvi a 9 ree each other. 
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The ability to construct a random interval from the value of 
o 

°W1 (and all other estimates) that contains the expected value 
of the backscatter with a known probability provides a new way to 
study the data set. If °vv^B’ and a vv^’ ^ do not 

in this random interval an appropriate fraction of the time, say, 
0.6826, if j 1 1 [ s |t 2 | = 1 is used, with due consideration of buoy 
contributions, then the model is not predicting the expected 
value of the backscatter and needs to be improved. 

Suppose, for example, that about 100 sets of data for re- 
ported buoy winds are collected. These would be 200 vertically 
polarized and 200 horizontally polarized sets of data for com- 
parison. The binomial distribution would be applicable. For two 
standard deviations, there ought to be somewhere between 124 and 
150 values from the model as calculated from and X g and, per- 
haps, from the MLE values 0 and X that fall within the ap- 
propriate random intervals. 

If much fewer than 124 successes occur, the result is strong 
evidence that the model has not correctly treated all of the 
various properties of the waves that produce backscatter and ac- 
counted for those effects that cause the backscatter to vary that 
are not directly related to the local wind speed and direction. 

If the number of successes is very much over 150, say, 175 
to 200, so that the random interval nearly always encloses the 
model values computed from U g and X g , then the random variation 
of the estimates of the backscatter is so great that there may be 
no skill in the recoveries of the winds from the backscatter 
data. The requirements for + 2 m/s or 10% for speed and corres- 
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ponding direction requirements will then have to be checked 
carefully. 

11.2 Comparisons of maximum likelihood estimates when con- 
ventional data are unavailable . Woiceshyn, et al. (1986) and 
Woiceshyn, et al. (1987) have uncovered many systematic sources 
of error in the wind recoveries from the SASS-1-S0S model and 
wind recovery algorithm simply from the study of the SASS wind 
data alone. Some way to compare various theoretical backscatter 
models without conventional wind data so that the best one could 
be further improved would be useful. 

Difficulties with the SASS-l-SOS wind recovery algorithm 
that were revealed by the study of Mode 4 data consisted of (1) 
systematic differences between the values of 0 for vertical 
polarization and horizontal polarization for high winds with the 
horizontally polarized wind speeds much higher (as much as 9 m/s) 
than the vertically polarized wind speeds and (2) systematic pat- 
terns in the wind direction recoveries such that for some wind 

speeds, wind directions for X(i.e. relative to x ) were dominant 

o 

for 45°, 135°, 225° and 315° and for other wind speeds only direc- 

. . o o o o 

tions of either 0 and 180 or 90 and 270 were dominant. The 

wind speed discrepancies need to be eliminated for a future model 

and the systematic direction biases would cause difficulties for 

any data assimilation method that depended heavily on correct 

wind directions. It would be possible to monitor the wind 

recoveries from a model and to determine areas of the ocean and 

wind speed ranges for which these various discrepancies occur. 
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Woiceshyn, et al. (1986) compared winds recovered from pairs 
of vertically polarized backscatter estimates 90° apart with 
winds recovered from pairs of horizontally polarized backscatter 
estimates 90° apart for the same location on the ocean for the 
available Mode 4 data from the SEASAT-SASS. The SASS-1 model func- 
tion was used with the SOS wind recovery algorithm, Jones, et al. 
(1982). Except for the logarithmic terms in (77), pairs of 
backscatter estimates about 90° apart often yield the same re- 
sults for both the SOS and the MLE. Differences are described in 
Pierson (1984) for the SASS-1 model function. Other slight dif- 
ferences will be investigated in terms of the design features of 
a scatterometer in a subsequent section. 

If the backscatter model that was chosen would have the 
property that horizontally polarized backscatter predicts higher 
winds than vertically polarized backscatter, the four backscatter 
estimates make it possible to detect this kind of discrepancy and 
to understand it so that the model can be corrected. Also consis- 
tently erroneous wind directions can be found and corrected. 

11.3 Detecting large wind speed discrepancies between 
models for vertically and horizontally polarized data . Although 
it is not essential to know the inverse of Eqn. (78) as in Eqn. 
(79) , enough can be learned about such inverses for any of the 
models described above to make it possible to describe the proper- 
ties of the backscatter estimates in terms of these inverses for 
analysis purposes. 
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Also if the sample values equal the model values, (77) is 
identically zero. Were it not for the terms involving the natural 
logarithms, the value of zero would be a maximum. The 
logarithmic terms are a logical consequence of the definition of 
a maximum likelihood estimate, and the actual maximum will be 
somewhat greater than zero and not at the value for an input 
without sampling variability. Maximum likelihood estimates are 
not necessarily unbiased. The effect of the logarithmic terms 
appears to be rather small, and it will be neglected in this sec- 
tion, to be discussed in more detail in terms of design data. 

Fig. 9a illustrates possible backscatter curves versus X for 
some fixed value of 0 for one of the outer 25 by 25 km areas of 
the swath. The figure is schematic so that no scale is given for 
the vertical axis. The upper continuous curve shows the variation 
of vertically polarized backscatter versus X for some fixed value 
Uj, for the incidence angle appropriate to the area. The in- 
cidence angles for the forward 45° beam and the aft 13^* beam are 
nearly equal so that the slight difference can be neglected for 
illustrative purposes. The next continuous curve from the top 
shows the variation of horizontally polarized backscatter versus 
X for Uj for the 115° beam for a lower incidence angle ap- 
propriate to that beam. The lowest continuous curve is for 
horizontal polarization versus x at the incidence angle for the 
45° beam for 0 .. The two dashed curves correspond to the values 
predicted by the model in use for U 1 for horizontal polarization. 

For the area under consideration, four estimates of the 
backscatter will be obtained. Suppose that the wind is from 
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Fig. 9 Model Errors that Result From an Incorrect Model for 
Horizontal Polarization: 

a) Backscatter Versus Aspect Angle for a Fixed U for 
Vertical and Horizontal Polarization with the 
Continuous Curves Representing a Correct Model and 
the Dashed Curves Representing an Incorrect Model 
for Horizontal Polarization; 

b) Results for the MLE if the Incorrect Dashed Curves 
Were Correct; 

c) Results for the MLE if the Correct Curves are the 
Continuous Curves for Horizontal Polarization and 
the dashed incorrect curves are used. (See Text) . 


c--* 
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o 

150 . The forward vertically polarized estimate will have an ex- 
pected value given by the point labeled, 1, and will depart in a 
random way from the point shown according to Eqn. (12) . The for- 
ward horizontally polarized value will scatter in the same way 
about the point labeled, 2. The point labeled, 3, is for the 
ll^ 3 beam for horizontal polarization, and corresponds to a wind 
direction relative to that beam from 80°. The point labeled, 4, 
corresponds to the 135° beam and falls at a wind direction of 60° 
relative to that beam. The two points shown by the x's are for 
the model predictions for horizontally polarized backscatter. 

Fig. 9b illustrates what numerical results could be obtained 
under the assumption that the two lower dashed curves for 
horizontal polarization and the upper continuous curve for ver- 
tical polarization were correct. The figure is schematic and no 
vertical scale is needed. Consider the evaluation of the terms 
involving a yyi Ec l n * ( 77 ) • the continuous curve with a minimum 
at 0°, a maximum at 90° and so on, is the graph of the inverse of 
Eqn. (78) as in Eqn. (79) as X is varied. 

If the wind direction was from the direction, X, relative to 
X Q , then the wind speed would have to be the value graphed so 
that the corresponding backscatter value from Fig. 9a would have 
been the measured value. If the wind direction in Fig. 9a is 

changed, the inverse curve will change by moving up or down rela- 
tive to u r Relative to X Q for inverse curves for the 45° beam, 
all inverse curves have a minimum at 0° , a maximum near 90° , a 
minimum at 180°, and a maximum near 270° with similar properties. 
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properly phase shifted, for all of the other inverse curves. For 
error free data, they must intersect at the input wind speed and 
direction for a model and its inverse. 

The model values always equal the value of (with its 
sampling variability) on this curve. Thus whatever the estimate 
of c vvi» there will be a locus of points in the U, x plane such 
that these two terms in (77) are zero. Were (77) evaluated at 
points in the U , x plane not on the curve the numerical values 
would rapidly become more and more negative. The smaller the stan- 
dard deviation of the estimate the more rapidly the ridge-like 
form following the curve becomes very negative. 

The continuous curve with a maximum near 0° and a minimum at 

o o 
90 corresponds to the value labeled cr VV4 for the upper curve of 

Fig. 9a. Were (77) evaluated for the two values of °yyi and 
o^ 4 four peaks would appear at the four points where the two 
curves cross. These four values of wind speed and corresponding 
directions, one in each quadrant, represent the usual four vector 
winds recovered by the operational mode of the SEASAT-SASS. Only 
those points where the curves cross are consistent with a pos- 
sible vector wind for the area sampled. 

The dashed curve with the minimum of the two minima at 290° 
(360°-70°) represents the inverse curve for the horizontally 
polarized measurement from the beam at 115° . The minima for this 
curve are different because in Fig. 9a upwind horizontal polariza- 
tion is shown to be stronger than downwind horizontal 

o o 
polarization. The shift away from 90 to the forward beam to 70 

away from the forward beam produces the result that although this 


93 



dashed curve passes through the point corresponding to Uj at 
150° , it is far from the speeds and directions for which the pair 

of curves for the vertically polarized values crossed for the 

o o 

values near 30 and 210 . Two of the peaks xn the value of (77) 
would be greatly reduced by this third contribution to the 
evaluation. 

The remaining dashed curve is for the horizontally polarized 
estimate, a° v2 . It again passes through the value ^ a t 150°. 
The four curves in the vicinity of a direction of 330° come close 
together, but do not all pass through a common point. If (77) 
has a low value at this point, it could be quite a bit higher 
than the value at 150° and Uj. contours for Eqn. (77) for an ap- 
propriate area of the 0, x plane could be revealing for a chosen 
model and for known design values of a, 3 and Y . Had the measure- 
ments had no sampling variability and had the model been correct, 
Eqn. (77) would have its maximum at the assumed input, Oj, 150°. 

Fig. 9b also illustrate an important feature of Eqn. (77) . 
Sampling variability effects have been neglected. Each of the 
estimates for the area under analysis is subject to the effect of 
sampling variability as in Eqn. (12) . The random variables, say, 
tj to t have an equal chance of being either positive or 
negative, and thus each of the four curves in Fig. 9b has an 
equal chance of shifting up or down. There are 16 possibilities 
as shown in Table 4. 

If all four curves shift upward, they could all still inter- 
sect at nearly the same wind direction, but at a higher value 
than Uj. if all four curves shift downward, the result would be 
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TABLE 4. COMBINATIONS SUCH THAT a wl , o HH2 , 0^3 AND a W4 

CAN BE GREATER THAN ( + ) OR LESS THAN (-) THE 
EXPECTED VALUE IN THE ORDER GIVEN. 
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a lower wind at nearly the same direction. 


Other kinds of shifts could also result in points in the 0 , 
X plane such that the four curves come close together and produce 
a large value for the maximum of the maxima of the MLE. Large 
values for the maximum of the maxima of the MLE do not always 
mean that the speed and direction recovered by the MLE are more 
accurate than smaller values for the maximum of the maxima. 

The shifts in the curves as a result of these random effects 
can at times produce an incorrect wind speed and direction. For 
Fig. 9b, the cluster of crossing curves that produce intersec- 
tions near 330 as a result of shifts because of the random 
variability could move more closely together and perhaps even 
merge into one point. 

Since the dashed curves for horizontal polarization in Fig. 
9a are assumed to be incorrect, the pattern shown in Fig. 9b will 
not occur. 

If, in fact, the variation of horizontally polarized 

backscatter is shown by the two continuous curves in Fig. 9a for 

a wind speed of U 1 and if the model is using the values for the 

dashed curves for a wind speed of 0^ the result would be very 

different from the curves shown in Fig. 9b as shown in Fig. 9c. 

The effect of an incorrect model would be to shift the two curves 

for the horizontally polarized value for the inverse curves from 
° , o . 

°HH2 and °hh 3 upward in the U, X plane by many m/s because the 

horizontally polarized backscatter values that were actually ob- 
tained would be higher than the model predicts for Oj . Whatever 
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the maximum of the maxima of Eqn. (77) would be, it would have to 
be very small compared to values that would be obtained for (77) 
from Fig. 9b. There is no point in the 0, X plane where all four 
curves come close to the value Uj and 150° . This would be an 
operational clue that something was not being accounted for cor- 
rectly by the model. 

Without any knowledge of the actual wind from a conventional 
source, the two backscatter estimates for horizontal polarization 
could be used to recover (for this example, others would differ) 
four vector winds, and the two vertically polarized estimates 
could be used to recover four more. The differences in the 
speeds, which will be large, for the four horizontally polarized 
winds and the four vertically polarized winds provide data to 
permit adjusting the model so that the recovered winds will be 
more nearly in agreement. 

It could also be possible that the curve for vertical 
polarization is too high and would need to be lowered. A few con- 
ventional winds included in a set of such calculations can 
quickly provide information on which way the model needs to be 
adjusted especially if the random intervals derived above are 
used to determine how far off the model is compared to the de- 
sired result. The sampling variability of any new design needs 
to be known so that statistics can be calculated to determine whe- 
ther the effect is real or the result of sampling variability. 

There are four possible ways to pair the four estimates that 
would be obtained such that the two estimates that are used are 
either 70° or 90° apart for antenna pattern 1 in Fig. 2. Each 
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one can be identified in Fig. 9c for the example under analysis. 
Each pair would permit the rapid recovery of four wind vectors 
(or two or three as special cases) with, usually, one per 
quadrant relative to x • As in Fig. 9c, there would be 12 plus 2 
(usually 16) quite different speeds and directions. For each 

pair, the MLE speeds and directions can be used to compute the 
backscatter from the model that ought to have been the expected 
value for the other two backscatter estimates that were not 
used. If the random intervals defined by (34) , ( 50 ) and ( 51 ) , do 
not enclose these model values for an appropriate percentage of 
data sets for a given incidence angle, the data alone are suf- 
ficient to reveal the inadequacy of the model. 

11.4 Detecting model errors that result from incorrect 
upwind-crosswind and downwind-crosswind differences . Fig. 10 il- 
lustrates the reason for the SEASAT-A SASS-1 model function to 
have the propensity to recover the two types of wind direction 
recoveries that were unrealistic as described by Woiceshyn, et 
al. (1987) . It extends the concept to winds that could be 
recovered by this new instrument if the model used had similar 
errors. Fig. 10a illustrates two possible models for the varia- 
tion of backscatter with aspect angle. For vertical polarization 

at a fixed wind speed for the upper two curves, one is continuous 

o 0 

and one is dashed. The curves can be used for the 45 and 135 
beam. The next two, one dashed and one continuous, are for the 
115° degree beam for horizontal polarization, and the lowest two 
are for the 45° beam for horizontal polarization. 
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Fig. 


10 Model Errors that Result From Incorrect Upwind-Crosswind 
and Downwind-Crosswind Differences: 

a) Backscatter Versus Aspect Angle for Two Possibilities 
and the Four Expected Values for Each, 

b) Continuous Curves Correct and Used to Recover Winds, 

c) Continuous Curves Correct, Model Uses Dashed Curves 
for Expected Values, 

d) Dashed Curves Correct, Model Uses Continuous Curves 
for Expected Values. 
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The figures are schematic. No vertical scales are needed. 
It is assumed that the wind is Uj from 170° and four expected 
values are shown for the continuous curves and four for the 
dashed curves. The estimates will scatter about the expected 
values as in (12) . 

Fig. 10b is the result of assuming that the continuous 
curves is Fig. 10a are correct and are used to recover winds from 
actual data (with sampling variability neglected) . Fig. 10c is 
the result of assuming that the continuous curves correspond to 
the real data, and the model is the dashed curves. Fig. lOd is 
the result of assuming that the dashed curves provide the actual 
data, and the continuous curves are the model. 

The difference between the various continuous and dashed 
curves for Fig. 10a is that the upwind-crosswind and downwind- 
crosswind differences in the backscatter are greater for the con- 
tinuous curves than for the dashed curves . A larger upwind- 
downwind difference is also shown for the horizontally polarized 
aspect angle variation. 

Fig. 10b is what would be expected in the absence of sam- 
pling variability if the model were correct. The continuous 
curves are for vertical polarization and the dashed ones are for 
horizontal polarization. The greater upwind-downwind difference 
for horizontal polarization, plus the 115° degree beam, helps to 
make the one point in the 0 , X plane where all four curves pass 
through 0 and 170° distinct from other possibilities. 

Fig. 10c is the result obtained if the dashed curves in Fig. 


100 


10a are for the model and the continuous curves are for the ac- 
tual conditions. The "true" backscatter values for both vertical 
polarization and horizontal polarization at the assumed direction 
of 170° are too high compared to the model values for the ver- 
tically polarized beam at 45° and the horizontally polarized beam 
at 45°. The vertically polarized value for the beam at 135° (at 
80° for X ) and the horizontally polarized beam at 115° (at 100° 
for x) are too low compared to the model. 

When the model is used to try to recover the wind speeds and 

o 

directions, the inverse curves for the 45 beam, both horizontal 
and vertical have to move up relative to 0^ and the inverse 
curves for the 115° (horizontal) and 135° beam (vertical) have to 
move down. As shown in Fig. 10c, the upper and lower pairs of 
curves might not intersect at all. A maximum likelihood estimate 
based on the vertically polarized data alone would yield very 
weak maxima for (77) as would the two curves for horizontal 
polarization. The value of the maxima from (77) for all four 
input backscatter values would also be low. 

There would most likely be two maxima for (77) one near 0° 
and the other near 180° for the conditions shown. A true wind 
near 90° would place the lower curves in Fig. 10c well above the 
upper curves and give maxima for (77) near 90° and 270°. 

Horizontal lines are shown on Fig. 10a from the maxima and 
minima of the dashed curves to the point of intersection with the 
continuous curves. For any set of data recovered by NSCAT such 
that the wind direction for, say , the 45° beam lies within the 
range of these horizontal lines, the backscatter values obtained 
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by the scatterometer will be higher (or lower) that than those 
for the model because the up wind-crosswind and down wind- 
crosswind differences are larger than the model predicts. For 
this example, backscatter values larger than the model predicts 
for the upper continuous and dashed curves would occur. The range 
is + 25° to 30° about the maxima and minima of the various pairs 

of curves, and unless the wind direction is somewhere near 45°, 
135°, 225° , or 310° plus or minus the range for which the actual 

backscatter exceed the maxima and minima of the model, the result 
will always be like Fig. 10c, except that the curves will move 
closer together. For + 25° about the maxima and minima, 200 
degrees of possible directions will all be recovered by (77) as 
pairs (or single values for one of two directions) at either 0° 
and 180° or 90° and 270°. Even if backscatter values are 
recovered within the range of possible values required by the 
model, the values will correspond to wind directions closer to 
0°, 90°, 180° and 270°, and all recovered winds directions would 
be closer to 0°, 90°, 180° and 270° than the correct directions. 
Those points near where the continuous and dashed curves cross 
would be the only wind direction that would be nearly correct. 

Fig. lOd is the reverse of Fig. 10c. The continuous curves 
in Fig. 10a are the model values and the dashed curves are the 
actual values that could be recovered for some fixed wind speed, 
0 , as a function of aspect angle. For an assumed wind direction 
of 170° , the inverse curves for the 45° beams must move down 
because the actual backscatter values are lower than those for 
the model. Those for the 115° and 135° beams must move up 
because the actual backscatter values are greater than those 
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for the model. Eqn. (77) might recover a wind speed and direc- 
tion somewhere near 115° and 250° . The vertically polarized pair 
of backscatter values alone would yield four winds with direc- 
tions near 65°, 120 °, 245° and 290°. The horizontally polarized 

pair would yield winds with directions near 80°, 115° , 255° and 
o 

305 as sketched. 

From Fig. 10a, it would never be possible to obtain 
backscatter estimates that correspond to the values required by 
the model for the ranges of wind directions shown by the 
horizontal dashed lines. The ranges of wind directions that 
would be obtained by (77) must be for the open spaces between the 
horizontal lines of Fig. 9a. 

Even if a single vector, U(?), X (?) is selected by (77), it 
will be far from the correct value for an extensive range of pos- 
sible correct wind direction. The wind directions will be concen- 
trated within a small range of directions centered near 45° , 
135°, 225° and 315° . The illustration consequently explains the 
results of Woiceshyn, et al. (1987) in terms of the SASS-1 model 
function and pairs of vertically polarized or horizontally 
polarized backscatter estimates 90° apart from the SEASAT-SASS. 
For the SASS-1 model function, the disparity because of the model 
was even more concentrated near 45° , 135° , 225° and 315° than 
Fig. 10 illustrates. 

11.5 Additional considerations . Figs. 9 and 10 illustrate 
another important feature for attempts to develop a good model 
for backscatter from waves and to relate the backscatter to the 
local wind speed and direction. Any model can be used to 
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simulate backscatter values and the effects of sampling 
variability can be simulated by Monte Carlo methods. Since the 
model is also used to recover the winds, results such as Fig. 9c, 
10c or lOd cannot be simulated unless deliberate kinds of dif- 
ferences between the model used to generated simulated 
backscatter estimates and the model used to recover the winds are 
introduced. If this were done, there is hardly any need to in- 
troduce sampling variability. Curves such as those in Figs. 9 
and 10 quickly reveal the inconsistencies that will be the re- 
sult. 


Unless the effects of sampling variability completely 
dominate differences between the model used and what is obtained 
by a scatterometer, a subset of the data for incidence angles 
nearest to 20°, 25°, 30°, 35° and so on for the 45° forward beam 
for every fifth point parallel to the subsatellite track could 
quickly reveal the kinds of inconsistencies considered for Figs. 
9 and 10. The random intervals that can be computed for each of 
the four estimates per 25 by 25 km area can also quickly reveal 
whether or not unique functions such as Eqns. (7) and (8) exist 
as functions of U and X with the alternative that a much more com- 
plicated model needs to be considered. 

The pattern shown by Fig. 10c can be quickly detected by a 
validation procedure when a scatterometer becomes operational. 
The recovered winds will have directions near 0° and 180° as il- 
lustrated (or 90° and 270° as a second possibility) . The model 
values for the backscatter for each beam can be found for these 
winds. The four sample values each have an appropriate random 
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interval. If these random intervals do not enclose the model 
values (they may be very, very far apart) , then the first type of 
error illustrated by Fig. 10a is evident. 

In contrast, the pattern shown by Fig. lOd is more difficult 
to detect. The sample values occur for the model curves, but not 
at a correct aspect angle. A test similar to the one for Fig. 
10c would not reveal any inconsistency. A running count of 
recovered wind directions for selected incidence angles and 
recovered wind speeds could eventually show the concentration of 

directions near 45°, and so on, and the absence of directions 

o o 

near 0 , 90 and so on. 

Figs. 9 and 10 also illustrate an advantage of the sampling 
pattern that was used as opposed to the alternative described in 
Part 1.2. For the alternative pattern there would be only one in- 
verse curve with maxima at 90° and 270° . Two others would have 

. . o o , 

minima near 110 and 290 , and the last one would have minima 

near 90° and 270° . Three of the inverse curves would track the 
same general region of the 0, x plane. Only one inverse curve 
could intersect these three curves at four possible wind direc- 
tions. If this one curve is high or low as a result of sampling 
variability, there will be a correlation between the errors in 
the recovered wind speeds and the recovered wind directions. 
Also the upwind-downwind difference for horizontally polarized 

backscatter for high winds and high incidence angles for the for- 

o 

ward 45 beam for the illustrated sampling pattern could at times 
eliminate a false solution by patterns such as those in Figs. 9 
and 10. 
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Results similar to the above analysis were obtained by 
Pierson and Salfi (1982) for the study of the NOSS. There are 
numerous examples of curves similar to Figs. 9b and 10b. The 
kinds of errors found by Woiceshyn, et al. (1986) and Woiceshyn, 
et al. (1987) were not known at that time. 

11.6 Additional applications of the random intervals . If 
the logarithmic terms in (77) are omitted, and if , t^, and 

represent the sampling variability of the four backscatter 
random variables as in (12) , the substitution of four values for 
(12) into (77) yields 

MLE (0, X) = - (tj + t 2 2 + t 2 3 + t*) (80) 

which suggests that the MLE varies something like a gamma (or Chi 
Square) PDF. This is not correct for reasons given in section 
11.3. The minimization of (77) juggles these four quantities 
relative to the ’'true" wind by allowing one term to become larger 
if other terms can be made even smaller. 

However, the calculation of the MLE involves the separate 
calculation of each of the first four terms in the sum. These 
are usually combined to give single values for the maxima of 
(77) . The numerical values of the first four terms in (77) con- 
tain information that is not in their sum especially if the 
values are computed in terms of a standard deviation in the 
denominator before squaring so as to obtain the signs of the t's. 

The values from the model can be computed from each of the 
MLE estimates of 0, X but upon substitution into (77) the same 
values for the four tJs would result. Thus if one of the values 
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of is outside of the range + 1.645, the 90% random interval 
for the corresponding sample value will not include the model 
value. The MLE's will not be sensitive to some of the combina- 
tions of sign shown in Table 4, but they can be sensitive to 
forced maxima such that the four inverse curves do not come close 
to intersecting at a point in the U, X plane. 

For each of the four values of the t^, a success can be 
defined as a values within + 1.645 and a failure as a value out- 
side of this interval. The probability of a success if the "true” 
wind were known for a "true" model would be 0.90 and the 
probability of a failure if the "true" wind were known for a 
"true" model would be 0.10. For four values of the the re- 
sult can be interpreted in terms of the binomial distribution as 
in Table 5. 

For, say, 10,000 successive recoveries of maximum likelihood 
estimates, the successes will be biased high compared to the 
"true" values. This would decrease the probabilities for two or 
fewer successes. If, then, many more than 486 recoveries have 
only 2 successes, or if many more than 36 recoveries have only 
one success, of if there are many more than one recovery without 
a success, for all of the maximum likelihood estimates, further 

study of the model is needed. Repeated values for the maximum of 

2 

the first four terms in the MLE less than -4,(1.645) = -10.8241 
virtually assure that the model is incorrect. Also for all of 
the maximum likelihood estimates for a given cell, a simple count 
of the successes and failures provides a ranking for the wind 
speeds and directions that may be helpful in removing ambiguous 
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TABLE 5 PROBABILITIES FOR A BINOMIAL PROBABILITY DENSITY 
FUNCTION FOR A SAMPLE OF FOUR SUCH THAT THE 
PROBABILITY OF SUCCESS IS 0.90. 


P (4 

successes) 



0.6561 

P ( 3 

successes, 

1 failure) 

= 

0.2916 

P (2 

successes. 

2 failures) 

= 

0.0486 

P(1 

success, 3 

failures) 

= 

0.0036 

P(4 

failures) 


= 

0.0001 



Sum 

= 

1.0000 
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solutions, especially if patterns emerge as a function of the 
quadrant for each wind. 

11.7 Applications to the SEASAT-SASS data . Except for two 
features, the Mode 4 backscatter data from the SEASAT-SASS are 
similar to the data that would be obtained by this new design. 
The SASS data can be co-located in sets of four backscatter 
estimates, two for vertical polarization and two for horizontal 
polarization, that are about 90° apart. The cell centers do not 
occur as close together as those for the new design. Except for 
the availability of the beam at 115°, which often eliminates two 
of the four ambiguities, the SASS data are similar to the new 
data if a , $, y and 6 are known. 

However, any candidate model for backscatter can be tested 
against the above analysis and the model deficiencies, if any, 
can be found. There are ways to obtain a candidate model based 
on the SASS data for Mode 4 that can be compared with improved 
versions of the models described above. 
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12 RESOLUTION AND WAVENUMBER ALIASING 


12.1 Twenty five kilometer resolution . For 25 kilometer 
resolution, the scatter in the sample values of the four 
backscatter estimates will be similar to those shown in Fig. 5 
for any wind direction. The maximum likelihood values for the 
wind speeds and directions would scatter approximately as the 
distance between the 5% and 95% curves. For winds below about 20 
m/s the errors in the wind recoveries compared to data buoy winds 
would not be too great, especially if refinements are introduced 


for 

low winds. Any mesoscale variability 

as described 

in 

Part 

8.5 

would 

be masked by the 

effects 

of variability 

in 

the 

backscatter 

estimates. There 

would be 

eight terms 

to 

be 


evaluated in Eqn. (77) except that the logarithmic terms will be 
shown to be negligible in Part 13, to follow, so that they can be 
omitted. 

12.2 Fifty kilometer resolution . Averages of the values 
for two adjacent areas parallel to the subsatellite track can be 
formed and grouped to represent an area 50 km by 50 km. The ef- 
fects of sampling variability would be reduced as described 
above. With the logarithmic terms omitted, there would be eight 
terms in the applicable version of (77) . How much this would 
reduce the scatter in the wind recoveries could be determined for 
any candidate model function. 

12.3 One hundred kilometer resolution . From Fig. 5, the 
recoveries for high winds will not be very good for twenty five 
kilometer resolution. They may still not be good enough for 
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fifty km resolution. There is the possibility that horizontally 
polarized backscatter will not behave in the same way as ver- 
tically polarized backscatter and will continue to increase on 
past those speeds in these figures for winds as high as to be an- 
ticipated for the most intense extratropical cyclones in the 
historical record. If horizontally polarized backscatter behaves 
differently, more accurate high winds might be recoverable. The 
winds over the Antarctic Ocean from August to October may be the 
highest to be expected. 

Improved computer based numerical weather predictions are one 
of the primary uses for the scatterometer winds. The present 
numerical models have yet to achieve 100 km resolution in the 
horizontal, and this, or even coarser resolution, is adequate for 
these models except near fronts. The high winds in an intense ex- 
tratropical cyclone can easily extend over areas of 1500 kilometers 
in diameter, and a 100 km resolution for these winds would be an 
important contribution. 

Sixteen 25 by 25 km areas grouped in sets of four see (Fig. l) 
parallel to the subsatellite track would be the starting point. 
The values parallel to the track would be averaged and the four 
slightly different values of a, $ and Y would be averaged and 
then be reduced by a factor of 4. The MLE would then be computed 
by using 16 terms (see Part 13). What would result for averages 
of four at upwind was discussed for Figs. 5a to 5d. These 
ranges of uncertainty as in Fig. 6 are probably better than the 
winds reported by ships of opportunity (Dischel and Pierson 
(1986)). If the model should verify approximately as shown, the 
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most extreme winds will be recovered with considerable accuracy. 

12.4 Wave number aliasing and filtering . At any nominal scat- 
terometer resolution, there will be contributions to the estimates 
of the vector wind from both the sampling variability of the 
backscatter estimates and from the mesoscale variability for each 
of the actual ocean areas sampled by a given antenna. These pertur- 
bations will alias into lower wavenumbers. The usual assimilation 
procedures for winds in numerical models will suppress these 
aliases to some extent. For some applications of these data, 
aliasing can be reduced by considering 125 km resolution involving 
twenty five 25 by 25 km areas with the recovered wind located at 
the center square. The 125 km resolution areas could be overlapped 
by 62.5 km to produce a smooth wind field. Again, special atten- 
tion is needed near fronts. 

The various tests described for the 25 by 25 km areas can all 
be extended in a straight forward way to these coarser resolution 
wind recoveries. The random intervals become much shorter, and the 
tests will reveal whether or no the coarser resolution results are 
good. Many of the hypotheses and tests described above can be 
modified to include the possibility that locally the values of t R 
are nearly constant. 
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13. FURTHER CONSIDERATION OF THE MAXIMUM LIKELIHOOD ESTIMATES. 


13.1 Discussion . The use of maximum likelihood estimators to 
obtain estimates of the unknown parameters from a random sample 
with an assumed probability density functions is a powerful me- 
thod. For moderate and large samples, it often yields an estimator 
similar to the method of moments except that it is biased. Eqn. 
(77) is an application of the maximum likelihood concept in a way 
that is quite different from the usual applications of the concept. 


There are other ways to estimate 0 and X, and the superiority 
of the MLE for the model that will eventually be used is still to 
be proven. Tables 3.1 to 3.10 gave values for the MLE to be com- 
pared with o° = a° as another choice for the estimate. The dif- 

o 


ferences between 


MLE 


and 


M 


are quite small. 


13.2 Maximum likelihood estimates for a wider range of o° ♦ 
Given a random value of o° there are two possible estimates of 
the model values. One is o° L£ and the other is o°. The quan- 
tity 0 ^y; was computed in 5 dB steps from zero to -40 dB and 
expressed in dB for both a single estimate and averages of 4 
estimates with a, 3 and y reduced by a factor of 4. The values, 
D and D are the numerical results for 


U = a MLE(dB) ' °M(dB) 


(81) 
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o 

for single values of o and 

D=o° . a 0 ( 82 ) 

MLE(dB) M(dB) 

for averages of four values. These differences (or effectively 
ratios) are given in Table 6 for the values of a , 6 and Y in 

Table 2. For high relative values of a 0 , the differences between 
the MLE and the expected value are small compared to the ranges 
of the random intervals that have been obtained. 

A strange result is that formally Eqn. (24) for low values 

o o 

of o can have a maximum for a negative value of o M . This can be 

shown by substituting a 0 = o into (24) . The natural logarithm 
becomes (1 + (3 /y)o° + (a/Y)(o°) 2 ) and for slightly negative 

values of the 3/y term dominates. The theory is then con- 

tradicted by the physics and the alternate choice, mentioned 

o^o 

above, namely ^ MLE - 0 = o M is correct. 

Eqn. (77) has the disconcerting property that one of its 

maxima will not be exactly at the input wind speed and direction 

if the random variables are replaced by their expected values 

throughout from a model. Since the values for the maximum 

likelihood estimates differ in almost a negligible way from the 

expected values and since the random interval for a given sample 

. o 

value are much much larger than the difference between 0 mle an< * 
it seems that the terms in Eqn. (77) involving the logarithms 
can be dropped in the search for the winds that fit the model. 

The new modified estimator would need a new name. What 
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TABLE 6a VALUES OF THE DIFFERENCES D = ° M LEfdB) " °MfdBl FOR THE RANDOM VARIABLE 
IN DB FOR INCIDENCE ANGLES, 0 IN DEGREES, FROM 15.6° TO 59.1°. 

NEG = MLE NEGATIVE. 
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should it be? Can it be proved to be superior to Eqn. (77)? Are 
the differences between the winds recovered by (77) and by (77) 
with the logarithmic terms omitted so small that the superiority 
of one over the other cannot be demonstrated by comparisons with 
conventionally measured winds? Answers to these questions are 
still to be determined. 

The important and basic feature of (77) is still retained if 
the logarithmic terms are omitted. The denominators of the four 
remaining terms, and as extended to lower spatial resolution all 
similar terms, are a function of the model values and not the 
sample values. This assigns the proper weights to each term in 
the sum and places greater emphasis on the terms with the lowest 
variance. 
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14 SUMMARY 


14 . 1 The probability density function for backscatter as a 
random variable . The PDF for backscatter estimates is very dif- 
ferent from the usual normal PDF in that the variance is a known 
function of the first moment (or mean) . This feature is used in 
several different ways to obtain statistics, with appropriate 
probabilities, for the interpretation of models for the expected 
value of the backscatter as a function of polarization, incidence 
angle, aspect angle and wind speed. 

The backscatter values obtained by a scatterometer are random 
variables. These random variables may either be completely indepen- 
dent in the probability sense from one estimate to the next or par- 
tially correlated because the unknown errors for the calculation of 
o° as a result of orbit and attitude uncertainties are slowly 
varying along the orbit. A test to find out whether or not succes- 
sive values of o° are truly independent, or partially correlated, 
has been devised. 

Statistical estimation procedures require that all calcula- 
tions be based on sample values. It is usually also required that 
the sample values be truly random and independent. Time series 
theory, which involves spectra and cross spectra, is needed if the 
successive random variables are not white noise; i.e. independent 
in the probability sense and normally distributed. 
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Quantities calculated from random variables that do not in- 
volve the unknown parameters of an assumed probability density 
function can be used to determine either point estimates for the 
unknown parameters or intervals with random end points with known 
probabilities of containing the unknown parameters. Such quan- 
tities are statistics. All of the above analysis is based on the 
calculation of various statistics. 

14.2 Recovering wind speed and direction estimates . Given 
a model for the expected values of the backscatter as a function 
of polarization, wind speed, wind direction and incidence angle, 
the constraints imposed by the model are imposed on a random 
sample of backscatter estimates so as to maximize the probability 
that the independent sample values came from PDF's with expected 
values and variances determined by the model. The result is 
similar to a maximum likelihood estimate in form, but it is based 
on products of at least four different PDF's. 

The resulting equation contains terms that involve the 
natural logarithm of the variance, which would follow from an 
exact application of probability theory. However, in general, 
maximum likelihood estimators can be biased. Calculations of the 
maximum likelihood estimates for one sample and for the average 
of four samples show that the maximum likelihood estimates differ 
negligibly from the expected value of the random variable. It is 
concluded that the logarithmic terms have a small effect on the 
estimates of 0 and X and that they introduce other undesirable 
theoretical features. They can be left out of the equation. 


119 


15 ACKNOWLEDGEMENTS 


This research was supported by the National Aeronautics and 
Space Administration (contract NAGW-690) and by the Jet Propul- 
sion Laboratory (contract 957714-02/CIT/JPL) . Dr Chong-Yung Chi 
provided the data for Tables 1 and 2. 


120 



16 REFERENCES 


Bahar, E. { Review of the full wave solutions for rough surface 
scattering and depolarization: Comparisons with geometric 

and physical optics, perturbation and two scale hybid 
solutions. J. Geophys. Res. 92, C5, 5209-5225, 1987. 

Banner, M. L. , and E. H. Fooks, On the microwave reflectivity 
of small-scale breaking water waves. Proc. R. Soc. London, 

Ser. A, 399, 93-109, 1985. 

Chi, C. Y. , D. G. Long and F. Li, Radar backscatter measurement 
accuracies using digital Doppler processors in spaceborne 
scatterometers, IEEE Trans. Geoscience and Remote Sensing 
GE-24 . 3, 426-437, 1986. 

Chi, C. Y. , D. G. Long and F. Li, Round-off noise analysis for 
digital signal power processors using Welch's power spectrum 
estimation. IEEE Trans on Acoustics, Speech and Signal 
Processing, ASSP-35 . 6, 784-795, 1987. 

Dischel, R. s., and W. J. Pierson, Comparisons of wind reports 
by ships and data buoys. Proc. MDS '86, Marine Data Systems 
International Symposium, Marine Tech. Soc. Gulf Coast Section, 
447-466, 1986. 

Donelan, M. A. , The dependence of the aerodynamic drag coefficient 
on wave parameters, in Proceedings of the First International 
Conference on Meteorology and Air Sea Interaction of the 
Coastal Zone . 381-387, American Meteorological Society, Boston , 
Mass., 1982. 

Donelan, M. A., and W. J. Pierson, Radar scattering and equilibrium 
ranges in wind-generated waves with application to scatter- 
ometry, J. Geophys. Res. 92, C5, 4971-5030, 1987. 

Douglas, B. C. , D. C. McAdoo and R. E. Cheney, Oceanographic and 
geophysical applications of satellite altimetry. Rev. Geophys. 

25, 5, 875-880, 1987. 

Duffy, D. G. , and R. Atlas, The impact of SEASAT-A scatterometer 
data on the numerical prediction of the QE-II storm, J. 

Geophys. Res., 91(C2), 2241-2248. 1986. 

Durden, S. L. , and J. E. Vesecky, A physical radar cross section 
model for a wind driven sea with swell, IEEE J. Oceanic Eng., 
OE-10 (4 ) , 445-451, 1985. 

Fischer, R. E. , Standard deviation of scatterometer measurements 
from space. IEEE Trans. Geosc. Electron, GE-10 , 106-113, 

1972. 


121 


Gilhousen, D. B. , An accuracy statement for meteorological 
measurements obtained from NDBC moored buoys. Proc. 

MDS '86 Marine Data Systems International Symposium, 

Marine Tech. Soc. Gulf Coast Section, 198-204, 1986. 

Guissard, A. C. , Response to "Review of the full wave solutions 
for rough surface scattering and depolarization: Comparisons 
with geometric and physical optics, perturbation and two- 
scale hybrid solutions. J. Geophys. Res. 92^, C5, 5225-5228, 
1987. 

Hasse, L. , On Charnock's relation for the roughness at sea in 
Monahan, E. C. and G. MacNiocaill (eds) Oceanic Whitecaps 
and their Role in Air-Sea Exchange Processes , 49-55, D. Reidel 
Publishing Co., 1986 . 

Iwata, N., An algorithm of microwave backscattering from a 
perturbed sea surface in Toba Y and H. Mitsuyasu, eds. , 

The Ocean Surface: Wave Breaking, Turbulent Mixing and 

Radio Probing ; 283 - 286 , D. Reidel Publishing Co., 1985 . 

Jones, W. L. , L. C. Schroeder, D. H. Boggs, E. M. Bracalente, 

R. A. Brown, G. J. Dome, W. J. Pierson, and F. J. Wentz, 

The SEASAT-A satellite scatterometer : The geophysical 
evaluation of remotely sensed wind vectors over the ocean. 

J. Geophys. Res., 87 (C5) , 3297-3317, 1982. 

Kwoh, D. S., and B. M. Lake, A deterministic, coherent, and 
dual-polarized laboratory study of microwave backscattering 
from water waves, I, Short gravity waves without wind, IEEE 
J. Oceanic Eng., OE-9(5), 291-308, 1984. 

Lyzenka, D. R., A. L. Maffett, and R. A. Shuchman, The contri- 
bution of wedge scattering to the radar cross section of 
the ocean surface, IEEE Trans. Geosci. Remote Sens., 

GE-21 (4) , 502-505, 1983. 

Mood, A. M. , F. A. Graybeal and D. C. Boes, Introduction to the 
theory of statistics, Third Edition, McGraw Hill, New York, 

564 pp. 1974. 

Moore, R. K. , and F. T. Ulaby, The radar-radiometer, Proc. IEEE 
57, 4, 587-590, 1969. 

Pierson, W. J. and R. E. Salfi, Monte Carlo studies of ocean 
wind vector measurements by SCATT: Objective criteria and 
maximum likelihood estimates for removal of aliases, and 
effects of cell size on accuracy of vector winds, NASA 
Contract. Rep., CR-165837-1, 282 pp. 1982. 

Pierson, W. J., The measurement of the synoptic scale wind 
over the ocean, J. Geophys. Res., 88 (C3), 1683-1708, 1983a. 


122 



Pierson, W. J. Highlights of the SEASAT-SASS program: A review 

in Allan, T. D. ed Satellite Microwave Remote Sensing . 

Ellis Harwood LTD., Chichester England, 69-86, 1983b. 

Pierson, W. J., A Monte Carlo comparison of the recovery of 
winds near upwind and downwind from the SASS-1 model 
function by means of the sum of squares algorithm and a 
maximum likelihood estimator, NASA Contract. Rep., CR-3839, 

62 pp. , 1984 . 

Pierson, W. J., W. B. Sylvester, and R. E. Salfi, Synoptic scale 
wind field properties from the SEASAT-SASS, NASA Contract. 

Rep., CR-3810 , 152 pp. , 1984. 

Pierson, W. J., W. B. Sylvester and M. A. Donelan, Aspects of the 
determination of wind by means of scatterometry and of the 
utilization of vector winds data for meteorological forecasts, 

J. Geophys. Res., 91(C2), 2263-2272, 1986. 

Plant, W. J., A two-scale model of short wind-generated waves and 
scatterometry, J. Geophys. Res., 91 (C9), 10735-10749, 1986. 

Schroeder, L. C. , D. H. Boggs, G. J. Dome, I. M. Halberstam, W. 

L. Jones, W. J. Pierson, and F. J. Wentz, The relationship 
between wind vector and normalized radar cross section used 
to derive SEASAT-A satellite scatterometer winds, J. Geophys. 
Res., 87, 3318-3336. 1982. 

Valenzuela, G. R. , Theories for the interacation of electromagnetic 
and oceanic waves-A review, Boundary Layer Meteorol. 13, 61-85, 
1978. 

Wentz, F. J., S. Peteherych, and L. A. Thomas, A model function for 
ocean radar cross section at 14.6 GHz, J. Geophys. Res., 

89 ( C3 ) , 3689-3704, 1984. 

Wentz, F. J., L. A. Mattox, and S. Peteherych, New algorithms for 
microwave measurements of ocean winds: Applications to SEASAT 
and the special sensor microwave imager, J. Geophys. Res. 

91 (C3 ) , 2289-2309, 1986. 

Woiceshyn, P. M. , S. Peteherych, M. G. Wurtele, D. H. Boggs, M. 

M. Borowski, G. F. Cunningham, A. Davies, and G. Muttitt, 

Cover photo: SEASAT scatterometer marine wind analysis, 

IEEE J. Oceanic Eng., OE-10(4), 342 pp, 1985. 

Woiceshyn, P. M. , M. G. Wurtele, D. H. Boggs, L. F. McGoldrick, 
and S. Peteherych, The necessity for a new parameterization 
of an empirical model for wind/ocean scatterometry, J. 

Geophys. Res., 91(C2), 2273-2288, 1986. 


123 


Woiceshyn, P. M. , M. Wurtele, D. Anderson, D 
and P. Janssen, Wind field inference from 
backscatter: The problem as a whole IGARS 

87CH24329 Vol. II, 981-986, 1987. 


Boggs, G. Cunningham, 
ocean microwave 
'87, IEEE 


124 



NASA 

Natonai AefonauJcs and 
Space /Vyrnnsiraton 


1. Report No. 

NASA CR-4228 


Report Documentation Page 


2. Government Accession No. 


4. Title and Subtitle 

Probabilities and Statistics for Backscatter 
Estimates Obtained by a Scatterometer With 
Applications to New Scatterometer Design Data 


3. Recipient's Catalog No. 


5. Report Date 

April 1989 

6. Performing Organization Code 


7. Author(s) 

Willard J. Pierson, Jr. 


9. Performing Organization Name and Address 

CUNY Institute of Marine and Atmospheric Sciences 
The City College of New York 
Convent Ave. at 138th Street 
New York, NY 10031 

12. Sponsoring Agency Name and Address 

Oceanic Processes Branch 

National Aeronautics and Space Administration 
Washington, DC 20546-0001 

15. Supplementary Notes 


8. Performing Organization Report No. 

10. Work Unit No. 

11. Contract or Grant No. 

NAGW-690 

13. Type of Report and Period Covered 

Contractor Report 

14. Sponsoring Agency Code 
EEC 


16. Abstract 


The values of the NRCS, <r , obtained by a scatterometer are random variables 
whose variance is a known function of the expected value. The probability 
density function can be obtained from the normal distribution. Models for the 
expected value obtain it as a function of the properties of the waves on the 
ocean and the winds that generated the waves. Point estimates of the expected 
value, <r° are found from various statistics given the parameters that define 
the probability density function for each value. Random intervals are derived 
with a preassigned probability of containing the value, cr^. A statistical test 
to determine whether or not successive values of are truly independent is 
derived. The maximum likelihood estimates for wind speed and direction are 
found, given a model for backscatter as a function of the properties of the 
waves on the ocean. These estimates are biased as a result of the terms in the 
equation that involve natural logarithms, and calculations of the point 
estimates of the maximum likelihood values are used to show that the contribu- 
tions of the logarithmic terms are negligible and that the terms can be 
omitted. 


17. Key Words (Suggested by Author(s)) 

Scatterometer 
Air-Sea Interaction 
Radar Backscatter 
Surface Winds 

19. Security Classif. (of this report) 
Unclassified 


NASA FORM 1626 OCT 86 


18. Distribution Statement 

Unclassified - Unlimited 


Subject Category 32 


20. Security Classif. (of this page) 

21. No. of pages 

Unclassified 

136 


NASA-Langley, 1989 





